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5 : Abstract 
_^ 

We present results for one-loop matching coefficients between continuum four- 
fermion operators, defined in the Naive Dimensional Regularization scheme, and 
' staggered fermion operators of various types. We calculate diagrams involving 

gluon exchange between quark lines, and "penguin" diagrams containing quark 
^ ! loops. For the former we use Landau gauge operators, with and without 0{a) 

^ \ improvement, and including the tadpole improvement suggested by Lepage and 

Mackenzie. For the latter we use gauge-invariant operators. Combined with ex- 
isting results for two-loop anomalous dimension matrices and one-loop matching 
coefficients, our results allow a lattice calculation of the amplitudes for KK mix- 
ing and K ^ TTTf decays with all corrections of 0{g^) included. We also discuss 
the mixing of AS* = 1 operators with lower dimension operators, and show that, 
with staggered fermions, only a single lower dimension operator need be removed 
by non-perturbative subtraction. 
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1 INTRODUCTION 



Many of the low-energy effects of the electroweak interactions can be expressed as matrix 
elements of four-fermion operators between hadronic states. For example, CP-violation in 
K — K mixing is proportional to such a matrix element (parameterized by Bk)- Knowledge 
of these matrix elements is required to extract the parameters of the electroweak theory. 
Since the matrix elements involve non-perturbative physics, lattice QCD is well suited for 
their calculation.^ A necessary ingredient, however, is the relationship between continuum 
and lattice four-fermion operators. This "matching" can be done using perturbation theory, 
since the operators differ only for momenta close to the lattice cut-off. 

The general form of the relationship between lattice and continuum four-fermion opera- 
tors is 

= O't + E (7i? H^IH + c.,) of + 0{g') + 0(a) , (1.1) 

i 

where /i is the renormalization scale used to define the continuum operators, and a is the 
lattice spacing. In other words, to calculate the matrix elements of a continuum operator 
Qcont ^ we must use the linear combination of lattice operators given by the right-hand-side 
of this equation. 7^^'* is the one-loop anomalous dimension matrix (using the conventions 
of Ref. 10]), and is the finite part of the one-loop matching coefficient. The finite part 
depends on the continuum regularization scheme, on the type of lattice fermions used, and 
on the precise definition of the lattice operators. 

In this paper we calculate Cij for two classes of lattice four-fermion operators composed of 
staggered fermions. In general, such operators involve quark and antiquark fields at different 
positions, and must be made gauge-invariant. Our first calculation concerns operators which 
are defined by fixing to Landau gauge and dropping all gauge hnk factors. We compute the 
one-loop matching coefficients for such "Landau-gauge" operators coming from diagrams in 
which a gluon is exchanged between quark lines (see Fig. |]), which we refer to as X diagrams. 
These are the only contributions for operators transforming as 27's under flavor SU(3), e.g. 
the AS = 2 operator leading to K — K mixing, and the A/ = 3/2 part of the K ^ 
amplitude. We consider both the simplest Landau-gauge operators and "smeared" operators 
having no 0{a) corrections at tree- level. These operators have been used extensively in 
numerical simulations ||^, ^, ^ ^ , and the results have small enough statistical errors that it 
is important to include the perturbative corrections. 

Our second calculation is of the contribution of "penguin" diagrams to the matching 
coefficients (see Fig. At one-loop, the contributions of X diagrams and penguin diagrams 
add, Cij = cfj + cfj", and can be treated separately. Penguin diagrams are present for 
operators transforming as octets under flavor SU(3). Such operators contribute to the A/ = 
1/2 part of the K ^ utt amplitude, and in particular to CP- violation in these decays. 
Numerical calculations of their matrix elements are possible , although less advanced than 
those for 27's. A technical problem with the calculations is that there is mixing not only with 
four-fermion operators but also with bilinears of lower dimension. These must be subtracted 
non-perturbatively, e.g. using the method of Ref. P]. If one uses operators made gauge- 
invariant by the inclusion of gauge links ("gauge- invariant" operators), then the mixing is 

■^For a recent review using lattice matrix elements to constrain standard model parameters see Ref. W. 
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constrained by gauge invariance. In contrast, Landau-gauge four-fermion operators can mix 
with gauge non-invariant bilinears, and require many more subtractions. For this reason, we 
present resuhs for penguin diagrams only for gauge-invariant operators. 

In the continuum, the mixing of octet operators with bilinears of lower dimension is 
restricted by chiral symmetry and CP. It turns out that only a single such bilinear can 
appear, with just the overall factor undetermined. This result is crucial for the subtraction 
method of Ref. |^ to work. Complications arise with lattice operators, however, because 
of the fermion doubling problem. The issue is well understood for Wilson fermions, where 
the breaking of chiral symmetry allows mixing with a number of lower dimension bilinears 



[^, |T0|, |TT|, necessitating the use of more elaborate subtraction schemes. With staggered 
fermions, on the other hand, mixing is restricted by the axial U{1) symmetry. While a large 
number of lower dimension bilinears can appear, only one of these contributes to physical 
matrix elements. This allows one to use the subtraction method of Ref. The fact that 
only one operator need be subtracted has been either implicitly or explicitly assumed in 



previous discussions [0, 0, |Tj], but is demonstrated here. 

The mixing coefficients Cij depend on the choice of scheme used to define the continuum 
operators. In physical matrix elements this dependence cancels with that of the two-loop 
anomalous dimension matrix. Thus for our results to be useful they must be given for a 
continuum scheme in which the two-loop anomalous dimension has been calculated. Results 
for both 27 and 8 operators are available in both Naive Dimensional Regularization (NDR) 
and 't Hooft-Veltman (HV — Ref. [|15|) schemes. We choose to give results for the NDR 
scheme. 

This work extends previous studies using staggered fermions. Gauge-invariant four- 
fermion operators were discussed by Daniel and Sheard |]TB[, and the complete results for 



some of these operators were given by Sheard . The calculation for Landau-gauge four- 
fermion operators turns out to be simpler than that for gauge-invariant operators, requiring 
only the matching coefficients for bilinears. These have been given in our previous work 



henceforth referred to as PS. Ishizuka and Shizawa have independently calculated the correc- 
tions from X diagrams for unsmeared Landau-gauge four-fermion operators [^, , and our 
results are in complete agreement with theirs. This is a non-trivial check since the calcula- 
tions are done using somewhat different methods. They have also calculated the corrections 
from X diagrams for gauge-invariant four-fermion operators. No previous calculations have 
studied penguin diagrams. The complete one- loop results for gauge-invariant SU(3) octets 
can thus be obtained by combining from Ref. ||19[ with those for cfj" given here. 

A single staggered fermion corresponds to Nf = 4 Dirac fermions in the continuum limit. 



We follow Ref. |T^ and introduce a staggered fermion "species" for each continuum flavor, 
i.e. Xu, Xd, Xs etc. Thus the continuum theory to which we match has Nf degenerate copies 
of each quark in QCD. It is straightforward to relate the matrix elements in this continuum 



theory to those in continuum QCD |12|, |T3|, [14 



The outline of this paper is as follows. In sect. ^ we explain our notation for four-fermion 
operators, and discuss their symmetries. Section^ gives the details of the computations of the 
X diagrams for Landau-gauge operators, and sect. ^ the corresponding results for operators 
of phenomenological interest. In sect. ^ we present the calculation of the penguin diagrams, 
and give results for the matching of the continuum operators which contribute to kaon 
decays. We close in sect, ^with a discussion of the mixing with lower dimension operators. 
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Two appendices collect results on Fierz transformations and the matching coefficients for 
bilinears. 

Preliminary accounts of this work have been given in Refs. 



141 and 20 . Some of the 



numerical results in these two works are wrong, and are corrected here. 



2 NOTATION AND DEFINITIONS 

In this section we set up the notation required to specify the four-fermion operators, and 
discuss their symmetries. While this section is self-contained, it draws heavily on the pre- 
sentation given in PS, which should be consulted for further details. 

2.1 Bilinears 

We use a gamma matrix basis for both spin and flavor matrices. To enumerate this basis 
we use "hypercube vectors", four-vectors whose components are or 1, which are combined 
using modulo-2 operations. A general gamma matrix is labeled by a hypercube vector 

is = 7f^7f^7f^7f* • (2.1) 

We use Euclidean space gamma matrices, which are hermitian, and satisfy {7/i,7i/} = 25^,^. 
An alternative basis is built out of the complex conjugate matrices = 7^, 

= ^['^M'^!' , (2.2) 

where F is another hypercube vector. To keep the notation clear, we always use ^ matrices 
for flavor, and 7 matrices for spin. Useful shorthand notations are 

S = ll^rs = l,, 5= (1,1,1,1)^75 = 75 , 5 = (0,0,0,0)^75 = / , (2.3) 

7575' = 755' , ^5^5' = ^55' • (2.4) 

We label quark fields in a continuum theory with four degenerate fermions using upper 
case letters, e.g. Qa,a- Here a is the spinor index, and a the flavor index, both running from 
1 to 4. For the moment we suppress color indices. A general bilinear is specified by a spin 
and a flavor matrix 

O'st =Qa,alfCF'Qp,t = Qa,ahs®^Fr'^'Qp,t = Q{lS ® ^f)Q ■ (2.5) 

In the second step we combine spin and flavor matrices into a single 16 x 16 matrix [is^^f)- 
We use a variety of lattice transcriptions of these continuum operators. All make use of 



a set of matrices which are related to (75 ® ^f) by a unitary transformation ||T6|, ^ 



(75®eF)^B = iTr[7i757B7^] (2.6) 
= E ilUrils ® ^Fr'^'ihB)^' ■ (2.7) 

a,a,l3,b 
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In the new basis we have traded the indices {a, a} for a hypercube vector A. Perturbative 
calculations in momentum space require another unitarily equivalent set of matrices 



{is ® e 



f)ab 



CD 



\A.C 



F)CD 4' 



\D.B 



All three sets of matrices satisfy the multiplication rule exemplified by 



J2 ®^f)ab ilS' ® ^F') BC = i^SS' ® CfF' 



'AC 



(2.8) 



(2.9) 



B 



To define lattice bilinears we divide the lattice into 2^ hypercubes in one of the 16 possible 
ways. Points on the original lattice are then specified by a vector y labeling the hypercubes 
(with all components even), and a hypercube vector C determining the position within the 
hypercube. The 16 components of the staggered fermion field x for ^ given y are collected 
in to a single hypercube field 

xiy)c = -My + c) ■ (2.10) 

In the continuum limit this becomes equal to Q, when expressed in the appropriate basis 



22 



c 



Q{y)c = C2^c)^'Qp, 



Thus the lattice operator 



C>sFiy) = E xiy)c iis^^F)cD xiy)D 

CD 



Xils ® ^f)x ■ 



(2.11) 



(2.12) 



has the same flavor, spin and normalization as Og'p*' = Qi'ys ^ ^f)Q in the continuum limit. 
An important property of a bilinear is its "distance" , 



J2 ~ 



(2.13) 



For the lattice bilinear of Eq. |2.12| this is the number of links between the quark and 
antiquark fleld. From this we construct the "distance-parity" of the bilinear, (— l)'^, which 
is +1 if the distance is even, —1 if it is odd. This parity is useful because of the U{1)a 
symmetry, valid when ma = 0, 



X exp 



ia(75 0^5) 



ia(75 ® 6) 



X , X X exp 

where a is the rotation angle. Under this transformation, even distance bilinears rotate 



(2.14) 



Xils ® ^f)x cos a xhs ® ^f)x + isma xhs5 ® ^F5)x 



(2.15) 



while odd distance bilinears are unchanged. It follows that even and odd distance bilinears 
cannot mix, and that the renormalization of even distance bilinears related by multiplication 
by (75 ® ^5) is identical |^. It has been noticed in Ref. fl^ and in PS that, at one-loop, 
the latter result holds also for odd distance bilinears, although this does not follow from 
U{1)a symmetry. We now understand this result. As long as we consider diagrams in which 
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the X and x fields are not contracted, then we can do separate U{1)a rotations on the two 
fields. In particular, if the rotation angles a for x and x are opposite, then odd distance 
bilinears rotate as in Eq. p.l5| , while even distance bilinears are unchanged. It follows that 
odd distance bilinears related by multiplication by (75 (S> ^5) will renormalize identically. 

As mentioned in the introduction, we consider both gauge-invariant and Landau-gauge 
operators. In the former, we insert between x and x a matrix given by the average of the 
products of gauge links along the shortest paths connecting the quark and antiquark. In the 
latter, we fix to Landau gauge, and then use Eq. p.l2| without any gauge links. Landau- 
gauge operators have the advantages of being simple to implement numerically, and of having 
Feynman rules with no coupling to gluons. Their disadvantage is that, nonperturbatively, 
there may be Gribov copies, in which case the definition of the operators is ambiguous. In 
perturbation theory there is, however, no ambiguity, and we do not consider this problem 
further here. 

As shown in PS, both gauge-invariant and Landau-gauge bilinears differ from the corre- 
sponding continuum bilinears at 0(a). For Landau-gauge bilinears, however, it is straight- 
forward to "improve" the operators so that they have only O(a^) corrections at tree level 
1^ . One simply replaces the quark field in Eq. |2.12| by 



xiy)A ^ i E + 2A[1 - 2A^])a = ^,Y.xiy + A + 2/i[l - 2A^]) , (2.16) 

and performs a similar replacement for the antiquark field. The resulting "smeared" op- 
erators are spread over a 4^ hypercube. Since the fields are shifted by an even number of 
lattice spacings, smearing does not affect the spin-flavor content or the distance-parity of 
the operator. 

We also consider a further improvement of the various operators. Lepage and Mackenzie 
have shown that a large part of the difference between lattice and continuum operators is 



due to fluctuations of the gauge links, and in particular due to tadpole diagrams p4[. These 
corrections can be partially summed using a rescaled gauge link U'^ = U^/uq in both the 
action and the operators. As shown in PS, the prescription for staggered fermion bilinears 
is to multiply them by Mq "^ in order to match more closely with the continuum operators. 
Here nu is the number of gauge links in the operator (rif/ = for Landau-gauge operators). 
We call operators including this factor "t adpole- improved" . There are various choices for 
Uq] for reasons given in PS we choose Uq = {^Tt{U)l), the normalized trace of the link in 
Landau gauge. 

2.2 Four- fermion operators 

The general four-fermion operator is specifled by two spins and two flavors 



Xhs' ® ^F')x Xils ® ^f)x = ils' ® ^F') (75 ® ^f) ■ (2.17) 

The second form is a useful abbreviation for the operator which we only use when there is no 
danger of confusion with the product of two single-barred matrices such as that in Eq. |2.9| . 
The two bilinears in Eq. p.l7| can be of any type (gauge-invariant, smeared Landau-gauge, 



etc.), although we only use operators in which both are of the same type. For a given type 



5 



of bilinear, there are 256^ four-fermion operators, so that the matching matrices are huge: 
65536 X 65536! Fortunately we need only calculate small blocks of this matrix. In particular, 
as explained below, matrix elements of phenomenological interest can be obtained using 
operators for which S = S' and F = F', which we refer to as "diagonal" operators, and the 
associated operators (75 (8> ^f) (tss ® ^fs)- 

For certain calculations we can restrict our attention further, and consider only the 35 
linear combinations of the diagonal operators which transform as the identity representation 
of the Euclidean lattice rotation group.^ Such "diagonal scalar" operators suffice, for exam- 
ple, in a calculation of 0, [l^ • There are few enough such operators that we can display 



some explicit results, rather than simply giving numerical answers. 

To enumerate the 35 diagonal scalar operators we use a compact notation. Since the spin 
and flavor of both bilinears are the same, we need only specify the spin and flavor of one 
bilinear. The first 25 operators are those in which the spin and fiavor independently form 
scalars. For the spin we use the basis 

S = I -I , V = J2jf,-j^ , T =J2 It^ly ■ 7i.7M > 

^ = H 7/^75 ■ 757m , ^ = 75 ■ 75 , 

where the dot separates the matrices in the two bilinears.^ Note that the spin matrix in the 
second bilinear is the hermitian conjugate of that in the first. We use exactly the same basis 
for fiavor except that 7 — >■ ^. The first 25 scalar operators are obtained by combining the 
five possible spins with the five fiavors. Our notation is exemplified by 



[A X P] = (7m5 ® ^5) (75m ® es) , (2.18) 



X T] = E E (7m ® U) (7m ® ipu) ■ (2.19) 

M ^'<P 

In the continuum limit, these operators become Lorentz scalars. 

In the remaining 10 operators the spin and fiavor are not independent. We denote these 
operators as 

K X ^m] = E (7m®U (7m ®U , (2-20) 

M 

with X A^], [A^ X and [A^ x A^ defined similarly; 



[^M X ^m] = E (7m ® (7m ® - (7m ® ip-^) (7m ® ^s.m) . (2-21) 



and the analogous [T^ x 



MJi 



[^M X = E (7m5 ® ^M-5) (75m ® 6.m) - (7m5 ® ^pu) ® 6m) , (2.22) 



fJ,,U 



^Among the 65536 operators of Eq. 2.17 there are many other operators transforming as the identity 
representation, but they are off-diagonal, i.e. S' S and/or F' ^ F ||l7|] . 

^We are using 5* to denote both a hypercube vector representing all possible spins, and a scalar operator. 
Where confusion might arise, we denote the scalar operator by /. 
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and the analogous [T^ x A^]; and finally 



(2.23) 



with [T+ X T_|_] defined similarly.]^ Linear combinations of these operators are written 



[{V + A),xT,] = [V,xT,] + [A,xT,] . 



(2.24) 



These 10 operators are lattice scalars, but do not become Lorentz scalars in the continuum 
limit. 

Each four-fermion operator comes in two "color-types". The color indices can be con- 
tracted either between bilinears or within each bilinear 



This notation indicates the number of color loops when each of the bilinears is contracted 
with an external color singlet meson: operators of type (I) give rise to a single color trace, 
while those of type (II) produce two traces. 

Lastly, we consider the implications of the U{1)a symmetry for the mixing of four-fermion 
operators. For diagrams in which the fields in the operator are not contracted with each 
other (i.e. X diagrams), we can do independent U{1)a rotations on each of the four fields. 
Several results follow from this. 

• The distance-parity of each bilinear cannot be changed by mixing. Diagonal operators 
can be classified as even or odd distance because both bilinears have the same distance- 
parity. Thus it follows that even and odd distance diagonal four-fermion operators 
cannot mix. 

• The mixing coefficients for the operators (75 ® ^f) {1s5 ® ^fs) can be obtained from 
those for (75 eg) ^p) (75 (g) C^p) by applying an appropriate axial rotation to the second 
bilinear. This is equivalent to multiplying the second bilinear by (75 ® ^5) from either 
the right or left. 

• The corrections are also unchanged if we multiply both bilinears by (75 (S> ^5) from, 
say, the left. We call this operation V. This invariance is derived by applying an 
appropriate axial rotation to both bilinears. 

These results are valid to all orders in perturbation theory, up to corrections of 0{ma). We 
have checked them explicitly in our one-loop calculations. They are not valid for penguin 
diagrams, as we discuss in more detail in sects. |^ and ^. 



Ol = (7s ® ^F)abbs ® ^f),, ; On = (75 ® ^F)aai^S ® ^f) 



(2.25) 



3 CALCULATING X DIAGRAMS 



^The definitions of [T± x T±] diSer from those used in Ref. |14 by a factor of two. 
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3.1 Overview of method 



In this section we describe the calculation of the contribution to one-loop matching coeffi- 
cients from X diagrams. We do this only for Landau-gauge operators. As mentioned in the 
introduction, we ffist match operators in QCD onto those in a continuum theory with Nf 
degenerate versions of each quark. We consider the following operators in the latter theory 

= S{js ® ^f)D 5(75 ® ^f)D . (3.1) 

This choice of flavors puts the operator into a 27 of flavor SU(3), and ensures that only 
X diagrams contribute to the matching calculation. This operator will match onto lattice 
operators having the same flavor 



Oj = Xshs' ® ^F')Xd Xsils" ® ^F")Xd ■ (3.2) 

Throughout this section, and the following, we assume that such a choice of flavors has been 
made, but do not show the flavor indices explicitly. 

We consider only diagonal operators in Eq. because the continuum operators that 
we are interested in are diagonal in spin, and because actual numerical calculations take 
matrix elements between states having the same flavor F. There is no technical difficulty in 
extending the calculations to off-diagonal operators. The diagonal continuum operators do, 
however, match onto lattice operators which are both diagonal and off-diagonal, as indicated 
in Eq. |3^ . In particular, at one-loop, operators appear having F' ^ F", although S' = S". 
The matrix elements of such operators vanish, however, between lattice states of the same 
flavor, up to corrections suppressed either by powers of a, or by additional factors of g^. 
Thus we need only consider matching with diagonal lattice operators. 

We are also interested in the matching matrices for odd parity operators of the form 

Siis®^F)DSi^S5®^F,)D . (3.3) 

As explained in the previous section, however, we can obtain the matching coefficients for 
these with no extra work. 

The relation between continuum and lattice operators, Eq. |1 . 1| , is derived by comparing 
the matrix elements of lattice and continuum operators between external states consisting of 
two quarks and two antiquarks. The coefficients '-f^j^ and Cij do not depend on the external 
momenta or quark masses except through terms suppressed by powers of the lattice spacing. 
Thus we set the external momenta and the quark masses to zero. Infrared divergences are 
regulated by adding a gluon mass, which can be set to zero in the final result. Ultraviolet 
divergences are removed in the continuum using a variety of renormalization schemes. 

The diagrams which contribute at one-loop are those of Fig. |I|, with the gluon propagator 
being in Landau gauge. Following Refs. |16| and PS, we call the ffist three Xa, Xb and Xc 



diagrams, and the remaining two Z and ZT diagrams. For tadpole-improved operators the 
contribution of ZT diagram cancels with the 0{g'^) term in the expansion of Mq^, so one 
can exclude these diagrams from the start. This is an example of how the tadpole-improved 
operators are similar to continuum operators, for the latter are corrected by all except the 
ZT diagrams. 
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(Xa) 



(Xb) 



(Xc) 




(Z) 



(ZT) 



Figure 1: Diagrams contributing to renormalization of four-fermion operators for Landau-gauge operators. 
The two boxes represent the two bihnears which make up the operators. Only half of the X diagrams are 
shown; in the others the gluon exchanges are between the remaining two fermion lines. Similarly there are 
additional Z and ZT diagrams. 

We first consider the one-loop lattice matrix element. We take the operator to be of color- 
type (II), and explain at the end how the results are modified for color- type (I) operators. 

The Xa diagrams combine with the self-energy corrections Z and ZT just as in the 
calculation for bilinears. The corrections can thus be obtained from the results for Landau- 
gauge bilinears. One must simply add the results for each of the two bilinears contained in 
the four-fermion operator. We discuss this in more detail in Appendix p. 

It turns out that one can also obtain the results for the Xb and Xc diagrams using only 
those for bilinears. This is done in a number of steps. First, we note that to bring the Xc 
diagrams into canonical form one of the quark propagators has to be conjugated, since the 
gluon connects two quarks (or antiquarks). The fermion-gluon vertex remains unaffected, 
but conjugation reverses the momentum flow along the fermion propagator and thus changes 
its sign when m = 0. Z and ZT diagrams, which are unaffected by conjugation, can now be 
added to the Xb diagrams, and subtracted from the Xc diagrams, such that each calculation 
is like that for a bilinear. Schematically 



Note that these manipulations are possible because the color factors are the same for the Xb 
and Xc diagrams. 

For the Xb diagrams, the next step is to do a simultaneous Fierz transformation on 
the spin and flavor indices. This is necessary because the gluon connects a quark and 
antiquark which are part of different bilinears in the original form of the operator. The Fierz 



Xb + Xc = Xb- Xc^ = {Xb + Z + ZT) -{Xc + Z + ZT)^ . 



(3.4) 
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transformations for spin (or flavor) alone are standard. We use the notation 

r' 

where F is a general Dirac matrix. Combining two such transformations, we find 

S',F' 

For the diagonal scalars, which Fierz transform into one another, we give the Fierz table in 
Appendix ^ In general we use a computer program to calculate the Fierz table. Having Fierz 
transformed, we can use the results for the bilinear corrections, and then Fierz transform 
back to the original basis. 

For the Xc diagrams, an additional step is required. First, one of the bilinears is conju- 
gated using 

{1s^^f)ab = (Is^^Uba ■ (3-7) 

One then does a Fierz transformation, calculates the correction, Fierz transforms back, 
and finally conjugates again. All these steps are straightforward, but, given the number of 
operators involved, they rapidly lead to an explosion of terms. In practice we have written a 
computer program to do the various steps, and have checked it by working out some examples 
by hand. 

The same method works for the four-fermion operators of color type (I). The only dif- 
ference is that the Z and ZT diagrams combine with the Xb diagrams to form a bilinear 
correction, while one must add and subtract Z and ZT diagrams to the Xa and Xc diagrams 
respectively. 

For continuum operators the steps in the calculation are the same, except for one subtlety 
which we return to shortly. The end result is that lattice and continuum corrections combine 
just as in the calculation for bilinears. Thus what we require from our previous work is the 
bilinear matching coefficients defined by 

j^cont ^ j^u^t ^ c'p_|!_[2(l _ as) ln(7r//ia) + j + off-diagonal . (3.8) 

Here B denotes a bilinear, Cp = 4/3 is the color factor, and as, which depends only on 
the spin of the bilinear, is (4,1,0,1,4) for spin tensors {S,V,T, A, P). The coefficients Cj 
depend upon the spin and flavor matrices in the bilinear, upon its type (e.g. smeared vs. 
unsmeared), and upon the continuum renormalization scheme. Numerical values are given 
in PS. 

We do not need to keep the off-diagonal components of the one-loop mixing matrix 
in Eq. These cause mixing between operators having the same spin but different 

flavor. For example, the continuum vector current (^(7^ ® Cu)Q (/^ 7^ J^) matches onto a 
linear combination of the lattice operators xilfi ® ^i^)x cind x(7m ® ^/i)x (110 sum on /i). 
Such matching coefficients imply that continuum four-fermion operators constructed of two 
bilinears with the same flavor and spin will match in part with lattice operators in which 
the spins of the two bilinears are the same but the flavors differ. As discussed above, the 
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contribution of such operators can be dropped in the matrix elements between states having 
the same flavor. We note, however, that in a two-loop perturbative calculation one must 
retain such off-diagonal mixing at intermediate stages. 

We now return to the subtlety concerning continuum regularization schemes. Fierz trans- 
formation and charge conjugation apply only if the gamma matrices are in 4-dimensions, as 
in Dimensional Reduction (DRED) and Pauli-Villars (PV) schemes. For schemes involving 
n-dimensional gamma matrices, such as NDR, the method described above will not work. 
In such schemes, one has to enlarge the basis of operators to include "evanescent" operators 
which vanish when restricted to 4-dimensions. There is considerable ambiguity in how one 
defines evanescent operators, and different choices lead to different results for the finite parts 
of diagrams, and thus for the coefficients q. The method described above amounts to a 
particular, and unconventional, choice of evanescent operators. We denote the result with 
this choice as being in the NDR' scheme. To relate the results to the conventional choice of 
evanescent operators of Refs. [|], ^ we must make a finite correction. We explain this 



further in subsect. 4.5 



In the "correction" step of the above calculations we need to convert the corrections for 
bilinears into those for four-fermion operators, keeping only the diagonal part. For a general 
diagonal four-fermion operator of the form of Eq. |2.17| , we simply sum the corrections to 
each of the constituent bilinears. In practice we implement this numerically. For the 35 
diagonal scalars, which are linear combinations of operators receiving different corrections, 
the results are worked out in Appendix |B[ 

3.2 Color factors 

Each gluon exchange gives rise to a color factor 

8 

X! '^abTcd = l^adScb - l^abScd ■ (3.9) 
0=1 

In general this mixes operators of the two color types. Since the color mixing does not 
depend on the spin or flavor of the operator, it can be factored out. To display the result 
we use a vector notation 

= ( ] , (3.10) 



i(ii) 



where i only labels the spin and flavor of the operator. The matching equation can then be 
written 

icont ilat / < > < > < >\ — ^lat 

=0. + 4-2T.[M% +M% C, +M% CA O , , (3.11) 



167r2 

3 



where A^^j are the contributions from the Xa,b,c diagrams, respectively, together with 
associated self-energy diagrams, with color factors (and (^^/IGtt^) removed. Ca,b,c are the 
corresponding color matrices 



Cb-\ \ n 1 > C!c-\\ g 1 • (3-12) 
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3.3 An example: \V x S] 



To explain the method in more detail we calculate the lattice operator which matches onto 
the continuum operator [l^ x 5*]. This is part of the matching calculation needed to extract 
Bk- For simplicity of presentation, we adopt the following notation for rows of the matching 
matrices 

M\[VxS])^irZMtp-\ . (3.13) 

V / i=[Vy.S] 

3.3.1 Xa diagrams 

For the Xa diagrams, we can take the result directly from appendix 

M\[V ^ S])=2cvs[V X S] . (3.14) 



3.3.2 Xb diagrams 

We must perform the operations Fierz-correct-Fierz. Using the appendices we find 

IQM\ [V X S]) = {ess - CSP + 3Cyv2 + CVVO - 5CVA2 - CVA4){ [S X V] + [P X A]) 

+ i~2csv + 2csA - 2cvs + 2cyp)( [SxV]- [Px A]) 

+ {csS - CsP - 'iCvV2 - CvVO + 3CyA2 + cvaa){ [S X A] + [P xV]) 

+ {2csv - 2csA - 2cvs + 2cyp)( [SxA]-[Px V]) 

+ {csS + 6C5T + csp + 'icvv2 + Cyvo + 3cva2 + cvaa){ X 5] + [A X P]) 
+ {Acsv + 4C5A + Cvs + 3CyTl + 3CyT3 + Cyp)( [V X S] - [A X P]) 

+ {ess - 2csT + csp) [{V + A)xT] 

+ {cvs - CVTI - CVT3 + cvp) [{V - A) X T] 

+ {css + QCST + Csp - 3CVV2 - CyvO - 3CVA2 - CvA4){ [V X P] + [A X S]) 
+ i-Acsv - 4C5A + Cvs + 3CyTl + 3CyT3 + Cyp)( [V X P] - [A X S]) 

+ (ess -Csp) [Tx{V + A)] 
+ {-2csv + 2csa) [Tx{V- A)] 

+ (cvvO - CVV2 + CvA4 - Cy A2) [{V + v4)^ X T^] 
+ (Cyyo - Cyy2 - CyAi + CVA2) [T^, X {V + A)^,] 
+ (2cyT3 - 2cvTl) [T^ X (\/ - A) ^] 

+ off-diagonal 

We have simplified the answer by setting /i = vr/a, so that the anomalous dimension factor 
is absent. It is simple to reinstate this factor, however, as it depends only on the spin of 
the bilinear being corrected (Eq. |37^ ), and thus can be determined by replacing csf with 
2(1 — as) ln(7r//ia) in the above equation. For the present example, this gives 

M\ [V X 5];anom. dim.) = -6 ln(7r//ia)( [V xS]+[Ax S]) . (3.15) 

Note that no even-distance operators appear in A^f,? as required by the axial symmetries. 
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3.3.3 Xc diagrams 

The first step is to conjugate, which for [V x S] has no effect. Thus the ensuing Fierz- 
correct-Fierz steps give the same result as for the Xb diagrams. Finally, we must conjugate 
again, and multiply the total result by —1. This gives 

-16M% [V X S]) = {ess - CSP + SCVV2 + Cyvo - 3cyA2 - CyA4)( [S X V] - [P X A]) 
+ {-2csv + ^CSA - ^Cvs + 2cyp)( [S X V] + [P X A]) 

- {css - Csp - ScvV2 - Cvvo + 3CyA2 + CyA4)( [S X A] - [P X V]) 

- (2csv - 2csA 2cvs + 2cvp)i [S x A] + [P x V]) 

+ (css + 6c5r + csp + Scvv2 + cyyo + 3cyA2 + cvAi)i [V X S] - [Ax P]) 

+ {4csv + 4CSA + CvS + 3CyTl + 3CyT3 + Cyp) ( [V X S] + [A X P]) 

- {css - 2csT + Csp) [(V - A) x T] 

- {cvS - CVTI - CVT3 + cvp) [{V + A) X T] 

+ {css + 6C5T + Csp - 3cyy2 - Cyyo - "iCy A2 - CvAi){ \V X P] - [A X S]) 
+ (-4csy - 4.CSA + CVS + 3cyTl + 3cyr3 + Cyp)( [V X P] + [A X S]) 

- {css - Csp) [Tx{V- A)] 

- {-2csv + 2c5a) [Tx{V + A)] 

- (cyyo - Cyy2 + CyA4 " Cy A2) [{V - A) X T^] 

- (cyyo - Cyy2 - CyA4 + CVA2) [Tf, X {V - A) ^] 

- (2cyr3 - 2cyTl) [T^ X {V + A)^] 

+ off-diagonal 

The anomalous dimension term can be reconstructed from this result in the same way as for 
Mb 

M% [V X 5];anom. dim.) = 61n(7r//xa)( [V x S] - [A x S]) . (3.16) 
Once again, only odd-parity operators appear. 

4 RESULTS FOR X DIAGRAMS 

4.1 General features 

Using the method explained in the previous section it is straightforward to obtain the 256 x 
256 block of the matching matrices yW"'**'^ connecting diagonal four-fermion operators to each 
other. Both the Fierz transformations and the (diagonal part of) the bihnear corrections 
can be represented as 256 x 256 matrices. To obtain numerical results, we have written a 
computer program to multiply these matrices appropriately and put in the color factors. We 
have done this both for the full matrix and for the 35 x 35 block connecting diagonal scalar 
operators to each other. We have checked these programs by comparing them to each other, 
and to various analytical results such as those of the previous section. 

We do not need the full matching matrix if we take matrix elements between states of 
definite external fiavor. For example, in the mixing of [V^ x 5*], we are only interested in 
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the matching coefficients for the operators with flavor S. Operators with other flavors will 
contribute only at 0{g^a) or 0{g^). Because odd-distance operators mix only with each 
other, it suffices to consider the mixing of [V x S] with [A x S*] . Thus we want only a 2 x 2 
block of the matrices A^"'*'^. Including color types, we then have a 4 x 4 matching matrix. 
We give results only for such sub-matrices. 

We quote results for two types of operator: unsmeared and smeared, both tadpole im- 
proved. The results for operators without tadpole improvement are the same except that 
diagonal matching matrix elements are shifted: Ca Cu — 24.4661. This simple change 
follows from the fact that tadpoles appear only on external lines. 

To present numerical results we must choose a continuum regularization scheme. We want 
to use NDR, but as explained above, the method of the previous section leads to results in 
a slightly different scheme, NDR'. Since we have to add a finite correction to convert from 
NDR' to NDR, we might as well use an intermediate scheme for which the methods of the 
previous section apply, i.e. a scheme involving 4-dimensional gamma matrices. We choose 
Dimensional Reduction ^7.1 using the "easy subtraction" scheme of Ref. [||] (DREZ). The 



prime indicates that, for four-fermion operators, our method differs slightly from that of Ref. 



2q] . Another virtue of using DREZ is that the numerical results for bilinears are given in 
PS using this scheme. This allows the diligent reader to check the results given here more 
easily. 

It is straightforward to convert from DREZ to other schemes. Changing schemes intro- 
duces a spin dependent, but flavor independent shift in the bilinear matching coefficients 

Ci(scheme 2) = Ci(scheme 1) -|- t5(scheme 2) — t5(scheme 1) . (4.1) 

For the schemes we consider tj = tp and ty = tA- The values of ts, for DREZ', dimensional 
reduction using the subtraction scheme of Ref. |^ (DRED), naive dimensional regulariza- 
tion with an anticommuting 75 (NDR'), and Pauli-Villars using the subtraction scheme of 
Ref. @ (PV) are 

( 0.5, 0, 0.5) (DREZ') 

( 0.5, 0.5, 0.5) (DRED) , . 

( -0.5, 0, 1.5) (NDR') ^ ' 

{ 0, 0, 0) (PV) 

for spin tensors (/, V, T), respectively. These results are taken from PS, except that we have 
added the result for the PV scheme. The shift in Cj introduces a change in ^W'^'^ which 
depends on the spins of the four-fermion operators but not on their flavor. Including color 
factors one finds the change in the matching matrices to be 

Cij(scheme 2) = Cij(scheme 1) + ^ dfj [^^(scheme 2) — ^^(scheme 1)] , (4.3) 

S=I,V,T 

where the three matrices d^'^'^ are independent of flavor. We quote results for these below, 
in addition to those for Cij. 



''For the PV scheme we use Mg = \xy/e. for the gluon regulator mass. 
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A change in renormalization scale yU is a special case of changing schemes. Using Eq. |3.8| 
it is simple to show that 

^if = -64 + 24 . (4.4) 

We now present results for various operators that are used in present calculations. Fol- 
lowing PS, to give some idea of the size of the corrections we set fi = tt/ a and use a boosted 
coupling g"^ = 1.8. With these choices, the anomalous dimension contribution to matching 
vanishes, and the matching matrix is to be multiplied by g'^/lGir'^ ~ 1/90. 



4.2 [V X S] and [A x S] 

This is the part of the matching matrix needed for a calculation of Bk using external pseudo- 
Goldstone pions. We are interested in the submatrix of Cij with indices 

i,j = i[Vx S]j, [V X S]n, [A X S]j, [A x S]n) ■ (4.5) 

Using symmetry operation V (defined in the previous section), this matrix is the same as 
that with indices 



[A X P]j, [A X P]n, [V X P]j, [V x P]jj) . 



(4.6) 



The minus sign arises because, under V, V ^ —A. (The other transformations are S ^ P 
and T ^ T.) The sign is, however, irrelevant here, since it is an overall sign for the operators 
and does not affect Cij. 



The matrices needed for changing schemes are 





/ 9 


-3 


7 




1 








6 


-2 


6 


7 


3 


9 


-3 




I 6 


-2 





/ 



/ 7 3 

16 

-7 -3 

V -6 2 



-7 -3 \ 

-6 2 
7 3 

16 / 



, (4.7) 



so that 7'^'''' = —6d^. For the finite part of the coefficients, the results for the (tadpole- 
improved) Landau-gauge operators areQ 



c(unsmeared) 



c(smeared) 



f 2.0450 -0.6817 6.8467 2.9343 

5.8686 -1.9562 

6.8467 2.9343 2.5726 -2.2643 

V 5.8686 -1.9562 -4.2203 / 

/ 10.0757 -4.6508 12.5283 5.3693 \ 

-3.8768 10.7385 -3.5795 

12.5283 5.3693 10.2457 -5.1609 

V 10.7385 -3.5795 -5.2371 / 



(4. 



(4.9) 



These corrections are quite small — multiplying by 1/90 one sees that they are 10% or less. 
Without tadpole improvement the corrections would be somewhat larger. 



^The accuracy of these and following results is « ±0.0001. 



15 



Ishizuka and Shizawa have reported results for a variety of unsmeared operators, both 
Landau-gauge and gauge- invariant [0. They quote results for a rectangular block of the 
mixing matrix for Landau-gauge operators: Cj^, with i given by Eq. and k running over 
all operators. This contains the square matrix Eq. |4.8| . Using the conversion formula 

- eg* (Ref. [19]) = + (In TT + , (4.10) 

and accounting for our opposite definition of A and T, we find that the results agree within 
the quoted errors. 



4.3 [S X S], [T X S] and [P x S] 

For calculations of e', matrix elements of the operators [S* x 5], [P x S], [P x P] and [S x P] 
are required. The matching coefficients for these operators also involve [T x S] and [T x P]. 
Thus we present results for the submatrices having indices 



t,j = {[Sx S]j, [S X S]ji, [T X S]j, [T X S]n, [P x S]j, [P x S]n) 
which are equal to the submatrices having indices 

iJ={[Px P]j, [P X P]n, [T X P]j, [T X P]n, [S x P]j, [S x P]jj) 
For the anomalous dimension parts we find 



d' 







1 


21 


7 


3 


9 




-3 \ 











64 


6 


-2 












1 




42 


18 


54 - 


-18 


42 




18 




24 




36 - 


-12 








36 




12 








9 


-3 


7 


3 


1 




21 






V 








6 


-2 







64 ) 




/ 




3 - 


-1 


- 


-3 


1 \ 

































1 

























2 





























-3 


1 





3 


-1 










\ 














) 










/ 


27 


-9 


-7 


-3 


27 




9 \ 












-6 


2 










1 




-42 


-18 


10 


18 


-42 


-1 


8 


24 




-36 


12 





64 


-36 


12 






27 


-9 


-7 


-3 


27 




9 




V 








-6 


2 







; 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



(4.15) 
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while the finite parts are 



1.3864 


12.5267 
10.7371 
0.75 


5.9359 


23.8899 
20.4771 
0.75 

V 



-20.1799 
-59.1533 
5.3686 
-3.5791 
-0.25 


2.4296 
1.3531 
10.2385 
-6.8257 
-0.25 




2.0878 
1.7895 
6.7342 

- 
2.0878 
1.7895 

3.9817 
3.4128 
17.1791 


3.9817 
3.4128 



0.8948 
-0.5965 
-5.5910 
-10.7371 - 

0.8948 
-0.5965 

1.7064 
-1.1376 
-10.3096 
-13.7499 
1.7064 
-1.1376 



0.75 


12.5267 
-10.7371 
-9.3253 


0.75 


23.8899 
20.4771 
-9.2877 




-0.25 \ 


5.3686 
-3.5791 
11.9549 
26.5395 y 

-0.25 \ 


10.2385 
-6.8257 
12.4850 
28.1675 / 



(4.16) 



(4.17) 



respectively for unsmeared and smeared operators. 

These corrections are larger than those for the vector and axial operators. For the 
unsmeared operators, the range is roughly —0.65 to 0.3, which is sufficiently large that it 
may not be possible to do reliable computations. Calculations with smeared operators should 
be more reliable, since the corrections range only from —0.15 to 0.3. 



4.4 Flavor ^35 and ^345 operators 

Numerical calculations have also been performed with lattice pions having fiavors other than 
^5. For example, Ref. uses fiavors ^35, ^45 and ^345. The continuum operators that are 
required, e.g. 

E (7^5 ® 65) (75m ® ^53) , (4.18) 

are diagonal but not scalar. Thus we must use the full basis of the 256 diagonal operators. 
We present results only for fiavors ^35 and 1^345, as one can extract the results for all other 
fiavors using the V operation and the lattice rotation symmetry. 

In this subsection we give results only for the Cij. The matrices, which are independent 
of fiavor, have already been given above. The only subtle point is that one must reexpress 
them in the different bases that we use here. 

We begin with the results for the vector and axial operators needed for Bk- For external 
fiavor ^35 states, the continuum operators we need are 



^ = E (7^5 ® 65) (75/. ® 63)(/,//), E (7m ® ^35) (7m ® ^53)(/,//)] • (4.19) 
M ' M 

to match onto which we need the lattice operators 



j = [E (7m5 ® 65) (75m ® ^53)(/,//), E (7m ® ^35) (7m ® ^53)(/,//), 
M^3 Mt^S 



(735 ® 65) (753 ® ^53)(/,//), (73 ® 65) (73 ® ^53)(/,//)] • 
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The matching matrix for the unsmeared operators is 



2.5817 


6.8467 


2.9343 


3.9134 


-6 


2867 


6.8467 


5.0667 


5.8686 


-1.9562 


0.0000 


-14 


9467 


5.8686 


2.9343 


2.6723 


-2.5636 


6.8467 


2 


9343 


2.6358 


1.9562 


0.0000 


-5.0184 


5.8686 


-1 


9562 


0.0000 



6.8467 
V 5.8686 

while that for the smeared operators is 



-1.9562 
-2.4541 



(4.20) 



/ 10.2345 
0.0000 
12.5283 
V 10.7385 



-5.1274 12.5283 

-5.1477 10.7385 

5.3693 10.3081 

-3.5795 0.0000 



5.3693 10.7478 

-3.5795 0.0000 

-5.3481 12.5283 

-5.7363 10.7385 



-6.6672 12.5283 

-9.2539 10.7385 

5.3693 10.2013 

-3.5795 0.0000 



5.3693 ^^ 
-3.5795 
-5.0277 
-4.8819 / 



For flavor ^345 pions, we need continuum operators 



^ = E (7m5 ® C345) (75m ® C543)(7,77), E (^M ® C345) (7m ® C543)(7,/7)] • 

and the lattice operators 

j = lYl (7m5'»645) (75m '»^543)(7,77), E ill^ ® CsAb) ® ^543) ^jji), 

M=3,4 M=3,4 



(4.21) 
(4.22) 



E (7m5®C345) (75m ®C543 )(/,//), E (7m ® ^345) (7m ® 6 



M=l,2 

The matching matrices are 



;543, 



M=l,2 



/ 3.1770 - 

0.0000 - 
6.8467 

V 5.8686 - 

/ 10.4955 
0.0000 
12.5283 

V 10.7385 



-4.0775 6.8467 

-9.0555 5.8686 

2.9343 2.6615 

-1.9562 0.0000 

-5.9104 12.5283 

-7.2357 10.7385 

5.3693 10.2191 

-3.5795 0.0000 



2.9343 2.6615 

-1.9562 0.0000 

-2.5312 6.8467 

-4.9320 5.8686 



-2.5312 6.8467 

-4.9320 5.8686 

2.9343 3.1770 

-1.9562 0.0000 



2.9343 \ 
-1.9562 
-4.0775 
-9.0555 / 



(4.23) 



5.3693 10.2191 

-3.5795 0.0000 

-5.0812 12.5283 

-5.0245 10.7385 



-5.0812 12.5283 

-5.0245 10.7385 

5.3693 10.4955 

-3.5795 0.0000 



5.3693 \ 
-3.5795 
-5.9104 
-7.2357 J 



for unsmeared and smeared operators respectively. 

For scalar, pseudoscalar and tensor operators of flavor ^35 



(4.24) 



^ - [(/ (8) C35) {I ® C53)(/,//), (75 <H) C35) (75 ® C53)(/,//), 

matching requires the following operators 

j = [(/ ® 65) {I ® 63) (/,//), (75 ® 65) (75 ® C53) (/,//), 



C35) (7:^M ® ^53^ 



(4.25) 



E (73m ® 65) (7m3 ® E (73m5 ® 65) (75m3 ® 63)(/,/7)] 

M7^3 M7^3 
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The results for unsmeared and smeared operators, respectively, are 



/ -8.9254 
0.0000 
0.7500 
0.0000 
12.5267 
\ 10.7371 



/ -9.1536 
0.0000 
0.7500 
0.0000 
23.8899 
20.4771 



10.7553 
23.3405 
-0.2500 
0.0000 
5.3686 
-3.5790 

12.0827 
27.0947 

-0.2500 - 
0.0000 
10.2385 

-6.8257 



0.7500 
0.0000 
-5.7497 
0.0000 
12.5267 
10.7371 

0.7500 - 
0.0000 
-8.1815 

0.0000 
23.8899 
20.4771 



-0.2500 
0.0000 
1.2283 
-2.0648 
5.3686 
-3.5790 

-0.2500 
0.0000 
9.1665 

19.3181 
10.2385 
-6.8257 



2.0878 
1.7895 
2.0878 
1.7895 
6.8395 
0.0000 

3.9817 
3.4128 
3.9817 

3.4128 
17.2133 
0.0000 



0.8948 
-0.5965 

0.8948 
-0.5965 
-5.9070 
-10.8816 

1.7064 
-1.1376 
1.7064 

-1.1376 
-10.4124 
-14.0240 



Finally, for operators with flavor ^345, it is convenient to first 
and pseudoscalars, i.e. 



2.0878 0.8948 \ 
1.7895 -0.5965 
2.0878 0.8948 
1.7895 -0.5965 
5.2835 -1.2389 
0.0000 1.5667 / 

(4.26) 

3.9817 1.7064 \ 

3.4128 -1.1376 

3.9817 1.7064 

3.4128 -1.1376 
17.1090 -10.0995 

0.0000 -13.1896 j 

(4.27) 

present results for scalars 



i = [(/ (8) 645) {I ® C543)(7,77), (75 ® 645) (75 ® C543)(/,//)] 

which requires the lattice operators 



(4.28) 



3 = [(/ ® 645) {I ® 643)(7,//), (75 ® C345) (75 ® C543)(/,//)] , ^ {liJiv ® C345) {luix ® C543)(7,//)] ■ 

IJ,<v 

(4.29) 

The results for unsmeared and smeared operators, respectively are 



(4.30) 





/ 


-8.1285 


8.3646 


0.7500 


-0.2500 


2.0878 





8948 






0.0000 


16.9653 


0.0000 


0.0000 


1.7895 


-0 


5965 






0.7500 


-0.2500 


-8.1285 


8.3646 


2.0878 





8948 




I 


0.0000 


0.0000 


0.0000 


16.9653 


1.7895 


-0 


5965 






-8.8985 


11.3174 


0.7500 


-0.2500 


3.9817 


1 


7064 






0.0000 


25.0539 


0.0000 


0.0000 


3.4128 


-1 


1376 






0.7500 


-0.2500 


-8.8985 


11.3174 


3.9817 


1 


7064 






0.0000 


0.0000 


0.0000 


25.0539 


3.4128 


-1 


1376 


For the continuum tensor operator 













(4.31) 



IJt<v 



C345) {luiJ. 

the following lattice operators are required 

3 -■ 



C543) 



(4.32) 



[(/ ® 645) {I ® 643) + (75 ® 645) (75 ® 643))(/,77), E (7; 



645) (lun ® ^543^ 



(734 «) 645) (743 «) C543)(7,77), (712 ® 645) (721 <8) C543)(7,//)] 
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Here means the sum excluding the terms with fi = 1, 

results are 



2 and yU = 3, z/ = 4. The 



12.5267 5.3686 6.8395 -5.9070 7.1507 -6.8404 5.2835 -1.2389 
10.7371 -3.5790 0.0000 -10.8816 0.0000 -13.3707 0.0000 1.5667 



, (4.33) 



23.8899 10.2385 17.1940 
20.4771 -6.8257 0.0000 



-10.3544 17.2254 -10.4486 16.8368 -9.2828 
-13.8693 0.0000 -14.1205 0.0000 -11.0117 



(4.34) 

for smeared and unsmeared operators respectively. 

In summary, we see that the corrections for operators with flavors other than ^5 are of 
moderate size for all spins, not exceeding 20% for unsmeared operators, and 30% for smeared 
operators. 



4.5 Converting results to NDR 

To convert our results for Cij from DREZ' to NDR, we proceed in two steps. First, we convert 
to NDR', using Eq. 



Ci, (NDR') = Q, (DREZ') + ^ - 4 



(4.35) 



Note that this conversion depends only on the spin, and not on the flavor, of the operators. 
In the second step we convert from NDR' to NDR. We do this in the standard way by com- 
paring the one-loop matrix elements in the two schemes. Since this is entirely a continuum 
calculation, it also does not depend on the flavor of the operators. We use the notation of 



Refs. m,M 



c,j(NDR) = Cij(NDR') + Ar^^ . (4.36) 

For the bilinears, i,j = [{V - A)i, {V - A)ji], i,j = [{V + A)i, {V + A)ij], and i,j = 
[{S — P)i, {S — P)ii], only the Xc diagrams contribute, and the result is the same for all 
three cases 

-( ' 

-3 1 • 



Ar. 



(4.37) 



This notation means that the result applies for any flavor, e.g. i,j = {[V x S]j, [V x S]ij). 
For i,j = [{S + P)i, {S + P)n, Tj, Tjj], both Xb and Xc diagrams contribute and we flnd 



Ar,- 





/ 9 -3 


14 


M 


1 





12 


-4 


6 


-28 -12 


-5 


-9 




^ -24 8 


-12 


4; 



(4.38) 



This last result deserves some explanation. Previous calculations using NDR have con- 
sidered only four-fermion operators with spins V±A and S — P. To extend these calculations 
to S + P and T we must decide how to project against evanescent operators. In the notation 
of Ref. [^, we use the projectors 



Pi = 3^(1+75) ® (I+75) , P2 = lis E 7^.7.(1+75) ®7.7/.(l+75) 



(4.39) 



^1/ 
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where the sums run over n = 4 — 2e dimensions. Given a positive parity operator O, whose 
composition in terms of 5* + P and T is unknown, 



O = a{S + P) + b{T) + evanescent operators 
but whose projections are known, we use 



a 
b 



j_ 

36 





24 + 28e 



1^ 



-6 



18e 
13e 



P2{0) 



(4.40) 



(4.41) 



We have dropped terms of O(e^). 

We close this section with the example of Bk- When calculated on the lattice using the 
pseudo-Goldstone kaons as external states, this is the ratio of the matrix element of 



Ob = [{V - A) X P]i + [{V - A) X P]ij 



(4.42) 



between a K and a to the square of the matrix element of the bilinear A/^ = (745 ® S,^) 
between a K and the vacuum ||, p^ . 
Using the above, we find 



^cont 



(NDR) 



1 + 



9' 



IGvr' 



-(41n(7r/a/i) — 



(4.43) 



where 

50 = 10.6702 [A X P]j + 1.6598 [A x P]n - 10.1426 [V x P]j - 7.4626 [V x P]ji (4.44) 

for unsmeared operators, and 

50 = 13.1912 [A X P]i + 10.3174 [A x P]// - 13.0211 [V x P]^ - 12.1878 [V x P]^j (4.45) 

for smeared operators. In these results, the effect of changing from DREZ' to NDR is the 
coefficient —4/3 multiplying Og*. This conversion factor agrees with that of Ref. |ll9 |. 

For the denominator in Bk, the matching relation is the same in NDR and DREZ 
schemes 



cont \ 2 



(A 



lat\2. 



1 + 



9' 



IGti' 



;2CfCvs) 



(4.46) 



with 2CfCvs = 0, —3.8769 for unsmeared and smeared operators respectively. Without 
tadpole improvement, each of the numerical coefficients in both numerator and denominator 
is augmented by —24.4661. 

We can combine the corrections to the numerator and denominator oi Bk, since the latter 
are an overall factor. This amounts to subtracting the coefficient of g"^ in the denominator 
from each of the numbers in the results for the numerator. Doing this, we see that the 
results with and without tadpole-improvement are the same. This is because, for Landau- 
gauge operators, the tadpole diagrams are only on the external legs, of which there are the 
same number for both numerator and denominator. 
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(a) (b) 




(c) (d) (e) 

Figure 2: Examples of diagrams leading to mixing of four-fermion operators with bilinears. Each box 
represents a quark bilinear. 

5 PENGUIN DIAGRAMS 

5.1 Overview 

We now turn to diagrams in which one of the quarks in the four-fermion operator is contracted 
with one of the antiquarks to form a closed loop. Examples are shown in Fig. They give 
rise to mixing of four-fermion operators with quark bilinears (possibly containing gluon fields 
and derivatives). The mixing coefficients are proportional to a'^ where d = 6 and d' > 3 
are the dimensions of the four-fermion operators and bilinears, respectively. The appropriate 
treatment depends on d': 

• For d' < 6, i.e. mixing with lower dimension operators, the coefficients are proportional 
to inverse powers of the lattice spacing. Such ultraviolet divergent coefficients cannot be 
calculated to the required accuracy using perturbation theory [§. One must, instead, 
use a non-perturbative method to subtract any lower dimension operators that are 
produced. In the following section we provide a catalogue of such operators, and 
show that only one of them contributes to the matrix elements that are used in actual 
numerical calculations. 

• For d' = Q, a reliable perturbative calculation is possible. We perform a one loop 
computation in this section. 

• For d' > Q, the calculation is sensitive to infrared physics and thus unreliable. The 
mixing vanishes in the continuum limit, however, and can be ignored. 
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The one- loop diagrams are those shown in Figs. |^(a-d), together with diagrams of the 
same form having additional gluons attached to the bilinears. These latter diagrams add no 
new features to the following analysis, and we will not discuss them further. Note that we 
have chosen to attach the quark loop to a single bilinear. All continuum operators can be 
transcribed onto the lattice in this way, and we use this method both because it is simpler, 
and because it has been used in all numerical simulations to date. Given this choice, it turns 
out that the diagrams of Figs. 0(a) and (b) lead only to mixing with bilinears having d' < 6. 
The only one-loop diagrams giving d' = 6 operators are those of Figs. ||(c) and (d), and 
these we calculate in this section. They also give rise to mixing with a. d' = 5 operator, as 
discussed below. 

In the continuum a considerable simplification occurs upon using the equations of motion. 
In the chiral limit, the d' = 6 bilinears either vanish, or can be converted into four-fermion 
operators. This use of the equations of motion is equivalent to attaching a quark line to the 
gluon, converting Figs. |^(c) and (d) into the "penguin" diagrams shown in Fig. |^. On the 
lattice, not all of the d' = 6 bilinears that occur at one-loop can be rewritten as four-fermion 
operators using the equations of motion. Nevertheless, it turns out that the operators of 
interest can be so rewritten, so we couch our discussion in terms of a calculation of lattice 
penguin diagrams. We first do the general calculation, and then describe some applications 



in sect. 5.5 



We concentrate on gauge-invariant operators in this and the following section. Such 
operators receive corrections from both Fig. |5(a), in which a gluon connects the closed 
quark loop with a spectator quark line, and from Fig. ^(b), in which the gluon connects 
directly to the bilinear. The latter diagram is absent in the continuum; on the lattice it 
serves to cancel mixing with lower dimension gauge non-invariant operators induced by Fig. 
|^(a). This cancellation does not occur for Landau-gauge operators, for which Fig. ^(b) 
is absent. This illustrates the problem with using Landau-gauge operators for diagrams 
involving closed quark loops: they mix with gauge non-invariant lower dimension bilinears. 
This requires additional subtractions, and makes it difficult to extract phenomenologically 
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interesting matrix elements. We explain this further in sect. 

At higher orders, diagrams of the type shown in Fig. ^(e) lead to mixing with bilinears 
both of lower dimension and with a large number of additional operators having d! = 6. 
Although we have not calculated the mixing coefficients, we discuss in sect. EM the most 



interesting of the operators that arise. Of course, there is also mixing with d' = 6 operators 
due to X diagrams, which we have considered above. 

Finally, we note that the results of this section are unaffected by tadpole improvement, 
because one-loop tadpole diagrams are already accounted for when calculating X-diagrams. 

5.2 Notation 

Penguin diagrams occur for operators which transform as octets under the continuum flavor 
SU(3) symmetry (or which contain a part which so transforms). Aside from overall factors. 



which we discuss in sect. 5.5, the calculation is the same for all octet operators, so we 



consider in detail only the operator 



Os'F',SF = XsilS' CF')Xd XuilS ® ^f)Xu , (5.1) 

where color indices are suppressed. The continuum flavors have been chosen so that there 
is only one possible contraction leading to a penguin diagram, that having a u quark in the 
loop. Note that we consider the most general form of the four-fermion operator, in which the 
spins and flavors of the two bilinears are different. This requires no extra work for penguin 
diagrams because the first bilinear is a spectator in the calculation. 
The continuum operator corresponding to Os'f',sf is 

0^°,%^SF = S{-fs' ® ^F')D U{-fs ® ^f)U . (5.2) 

By construction, the tree level matrix elements of these two operators, calculated respectively 
on the lattice and in the continuum, are the same for physical momenta^ 

{0's'f\sf)^'^ = {Os'F',SFy'\l + O(pphysa)) . (5.3) 

The superscript refers to the number of loops. Both lattice and continuum theories differ 
from QCD by having an additional factor of Nf = 4 flavors. In particular, there are Nf 
quarks running around the loop, rather than a single quark in QCD. To match onto QCD, 
one must divide the diagram by Nf. To make it possible to see the effect of this correction, 
we keep the factor of Nf explicit below. 

One-loop diagrams lead to mixing with a single class of SU(3) octet operators, which we 
call "penguin operators". The continuum and lattice forms are, respectively. 



'PS'F',S"F" = XsilS' '^^F')Xd J2XqilS" ®^F")Xq , 

q 

V^f,^s"F" = S{js'^^F')D Y.Q{^s"®^F")Q . (5.4) 



^As discussed in PS, this is true when the continuum matrix element is expressed in the appropriate 
spin-flavor basis. 
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The sums run over the number of active hght flavors, typically m, d, s and c in present lattice 
simulations. 

The color factors are simple to calculate. All the operators come in two color types, and 
the matching matrix is a 2 x 2 matrix in color space. Using the notation of sect. P]^, we find 

tcont >lat < > > 

Oi =Oi +T7r^zZPii Cpcn P J + X diagram contributions , (5.5) 
where i and j represent both the spin and flavor indices of the operators, and 

Cpen= I f g g j (5.6) 



is the color matrix. The first row follows from Eq. |3.9| , while the second row vanishes because 
the gauge generators are traceless. 

To calculate p we need the lattice and continuum one-loop matrix elements. These have 
the general form (in a notation based on that of Ref. p8| ) 

{Os'F',Sf)% = E PfF,S"F"{Vs'F',S"F")^'^ (5.7) 



S"F" 



(the subscript pen indicating that only the penguin part of the matrix element is included), 
and 

/(■ncont \{1) _ y cont /jycont \{0) /p- o\ 

\^S'F',SFlpen ~ ^ 2 2^ PSF,S"F"\' S'F',S"F"/ 5 
S" F" 

respectively. This is a rather elaborate notation for the continuum result, because, when 
ma = 0, almost all of the one-loop coefficients vanish 

P'sf!s"F" = E ^S^L ^FI ^S'^il ^F'l P^"j% ■ (5.9) 



Such a notation is required, however, for the lattice matrix elements. Combining Eqs. |5]7 

P^^=PT-P^ ■ (5-10) 



and pTq , and using Eq. we find 



This notation is slightly redundant in that pij depend only on S, F, S" and F", but not on 
S' and F'. 



5.3 Calculation 

Since the first bilinear i n Fig. ^ is simply a spectator, we are, in effect, calculating the mixing 
of the second bilinear, (75 ® C,f), into other bilinears. We first consider this mixing on the 
lattice. Using the Feynman rules given in PS, and the momentum definitions from Fig. |^, 
we find 

{hs^^F))% = E rei^s®^F)ABe-"'-^^^''^^' E ^.(^; A)D,^{k)V.{k) . (5.11) 



10 



In this subsection the sum over quark flavors q (see Eq. 5.4) for the lower quark line in Fig. ^ is implicit. 
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Here and in the following we have suppressed all color matrices and indices, since the color 
factors have already been calculated. The factors P, D and V represent respectively the 
penguin loop, the gluon propagator and the lower quark-antiquark-gluon vertex. Expressions 
for D and V are given below. We first focus on the penguin loop, for which we find 



P^{k-B,A) 
h 
h 



4 



/ (/i - l2)e 



i(t>-{A-B) 



CD 



J2^(^)bc c^(7m ® I)cd^(^)da 

CD 



(5.12) 
(5.13) 

(5.14) 



where is the quark propagator 
S{(I))ab 



AB 



m 



(5.15) 



and we are using the notations 



= 167r' 



U (27r)4 

Various features of these results require explanation. 



= cos( 



(5.16) 



The factor of 1/16 multiplying (75 ® ^f) in Eq- |5.11| arises from the fact that we 
integrate the momenta over the entire Brillouin zone and trace over the spin-flavor 
indices. To avoid "double" counting we have to divide the result by 16. 

Infrared regularization is provided by the quark mass, which we can only set to zero 
after the lattice and continuum results are combined. 

We have transformed the matrices to the single-bar representation, which corresponds 
to physical position on a hypercube. This is the most convenient prescription since 
it allows one to easily determine the constraints on A and B. We have used these 
constraints to simplify Eq. |5.12 . 

The I2 term in Eq. |5.12| is the contribution of Fig. |^(b). We have written it so as to 
show explicitly that vanishes when k = 0. 



This vanishing is one of the consequences of the Ward identity, 

Y,sm{lk^)P^{k;B,A)=0 , 



(5.17) 



which is satisfied by our expression for all values of k. This identity follows from gauge 
invariance, and is necessary for the gluon to remain massless. 
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It turns out that contributions to operator mixing come not only from A; ~ 0, as in the 
continuum, but also from ~ Ctt, where C is a non-zero hypercube vector. This is because 
a momentum close to tt is physical for staggered fermions. We first consider A; ~ 0. The 
integrand Ji has to be expanded out to 0{k'^). Tedious algebra leads to the result 



Pf,{k; B, A) = ie^p„^kp{-i„^ ® ^^5) - i\mkp (7^^ ® I)^^ - (7^^ (g) I) 



'BA 



I, 



b 



Nf 



Ir = ^ IF^ 



HKku - ® + 0(P) , (5.18) 

where the indices p, a and u are summed. The integrals are 

/ F^ldsl = ^ X 11.2291 , h = ^l Faciei = ^ {j - 40.7767) , 

2 - 2sl - sls"^ =^i^^j F^- 9.5146) . (5.19) 

Note that la is finite, while and Ic diverge logarithmically as m — > 0. 

The first two terms in Eq. [5.18| correspond to mixing of the bilinear with gluonic operators 
which cannot be rewritten using the equations of motion. In other words, for these terms we 
are calculating Figs. |^(c) and (d) rather than penguin diagrams. For these terms we should 
drop the factor of x y in Eq. ^.11| , and multiply by the gluon polarization tensor e^. 

To analyze the term we must expand the factor e-«'= (-4+s)/2 jj-^ ]^q_ ^.11| up to linear 
order in k around the center of the hypercube. The leading term is unity, and gives rise 
to mixing between the bilinear x{lau ® ia5)x and the lower dimension gluonic operator F^i, 
(the dual of the gluon field strength). Reinstating the spectator bilinear, this is an example 
of mixing with a ci' = 5 operator. It is consistent with the lattice symmetries, but, as 
explained in sect. |^, it does not contribute to the matrix elements of practical interest. For 
the 0{k) term in the expansion, one finds after some algebra, that it leads to mixing between 
x(7cr5 ® icTv)x a-nd DyFfj^. Including the spectator bilinear, this corresponds to mixing with 
a = 6 operator. At one-loop, this mixing does not affect the operators needed to study 
the AJ = 1/2 rule (which are listed below in sect. |5.5| ), since none contain bilinears having 
flavor = ^cru- It will affect these operators at two-loop order. However, as explained in 
sect. D, its contribution in matrix elements will be suppressed by powers of a due to its 
flavor. 

The lb term corresponds to mixing of the operator xilf^p ® I)x (/^ 7^ p) with mF^p, 
both of which have dimension 3. Such mixing also occurs in the continuum, and when one 
combines the lattice and continuum corrections the infrared divergence in If, cancels with 
that in the corresponding continuum integral. Since the operator vanishes in the chiral 
limit, it is dropped in renormalization group analyses 0, and we follow this practice here. 
In fact, for the operators needed to study the AJ = 1/2 rule, the initial operator contains 
no terms with tensor bilinears, so mixing with mFpp begins only at two loop order. Thus its 
contribution is suppressed both by the quark mass and by an additional power of g'^/lGn'^. 
The factor of m appears because of the chiral symmetry of staggered fermions, and is not 
present for Wilson fermions |]lT] . 



We are thus interested only in the term proportional to Jc, which corresponds to the 
mixing of (7,^ ® I) with Dp^Fp^^, and which can be rewritten using the equations of motion. 
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To extract the infrared divergence we use 

1 F2 = 16(- In(4m2a2) - 7^ + Foooo) + OHmaf) . (5.20) 

Numerical values of the constants are given, for example, in PS. We calculate the remaining 
finite piece of Ic numerically, yielding, in total 

QJ\T 

h = Hma/n) - 8.8945 + 0((ma)2) . (5.21) 



The error on the numerical constant is roughly 0.0001. 

To complete the amplitude, the gluon must be joined onto the external fermion line. The 
result in Feynman and Landau gauges is the same because of the Ward identity Eq. |5.17| , 
so we use the Feynman gauge propagator 



For ~ the vertex is (see PS) 



V;(p) = -tgSip' + k- g')(7M ® ^A'B' cos(g; - ^k,) . (5.23) 

The external momenta have been decomposed into physical and flavor-related parts: p = 
p' + A'tc and q = q' + B'tt, where p', q' G [—n/2, vr/2]. We assume that p' and q' = k + p' 
are both close to zero, so that the cosine can be set equal to unity. Using the equations of 
motion on the vertex (i.e. X]/^ sin(A;^/2)V^ = for massless fermions), we see that only the 
k"^ term in the coefficient of Ic in Eq. |5.18| survives. This cancels the denominator of the 
gluon propagator, so that there is no longer any dependence on k, and we have arrived at a 
local operator. We can now do the integral over k to remove the delta function in Eq. |5.23| , 



and sum over A and B in Eq. |5.11| . We can set the factor exp{—ik ■ {A + B)/2) to unity, 
since the corrections lead to operators with d' > Q. Thus the only non-zero matrix element 
from gluon momenta A; ^ is 



Since (7^ I)a'b' ^'^^^ level matrix element of the bilinear (7^ ® I), we can read off 

one contribution to p'"* 

pfir,i = -^^ ■ (5-25) 

Note that this contribution is diagonal in spin and flavor. 

The corresponding continuum calculation is straightforward (see Ref. for more de- 
tails), with the result 

= -4^,ir'i:5srM(i,®r). (5.27) 
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The continuum integral is 

ln(— ), (5.2^ 



3 >2 



in DREZ, DRED, NDR, HV and PV scheme^^ Thus the only non-zero elements of the 
continuum mixing matrix are 

= -^c"""' ■ (5.29) 
The infrared logarithms cancel when we combine lattice and continuum results 

Pm/,^/ = - = Hf,a/ir) - 8.8945 . (5.30) 
The additional contributions discussed below do not change this mixing coefficient. 

5.4 Contributions from high momentum gluons 

As mentioned earlier, this is not the whole story. Unlike in the continuum, gluons with 
momenta k ^ Cn {C ^ 0) also contribute on the lattice. This is because, while p' and q' 
must be close to zero, k = q' — p' + n^B' — A') need not be. These contributions result 
in a local operator because both the gluon propagator and the fermion loop are both short 
ranged. 

It is instructive to begin by considering the vertex connecting the gluon to the external 

line 

V,{Cn) = -zgS{p' - q') ^^^^~^^^'^^ W^A',B'+c ■ (5-31) 

Clearly we must choose 0^ = for the vertex to be non-vanishing. Notice the appearance 
of the subscript _B' + C in place of B'. It is simple to show that 



V^iCn) = -ig5{p' - q'y- Y^(7.c ® ^c)a'B' ' (5-32) 

where C^, =2 J2p^u ^p- This shows that the external quarks "see" this vertex as having a 
different spin and flavor. 

The expressions f).31\ and ^.32| differ from the simplifled Feynman rules listed in PS. The 



rules given there are appropriate for propagators and vertices involved in loop integrations, 
but are inapplicable if more than one of the legs of the vertex is an external line. This 
is because, to derive the rules, one has to use the freedom to redeflne loop variables. For 
tree level couplings no rearrangement is allowed and the correct vertex is as given above. 
Note that similar changes occur for matrix elements of operators when they are inserted at 
non-zero momenta. 

The penguin loop is straightforward to calculate using the fact that the integration vari- 
able can be shifted freely. Furthermore, since the Ward Identity Eq. p.l7| holds for all k, the 
longitudinal part of the gluon propagator does not contribute, so that D^^, oc 6^^^. Thus in 

^^In the PV scheme /i is the fermion regulator mass. 
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evaluating we can assume that = 0. Including the momentum factor from the bilinear 
vertex, the loop becomes 



+M7,c®ec)^^E((-lf''-l)- (5-33) 
p^/i 

The integrals needed are 

/. = ^ / F^c?^^ = ^ X 16.3107 , Ie = ^J F'cjslsl = ^ x 5.0816 , (5.34) 

both of which are infrared finite. The first term in Eq. |5.33| corresponds to mixing with a 
d' = 7 operator, because of the factor of m, and can be dropped. The remaining two terms 
correspond to mixing with d' = 6 operators. 

To complete the calculation we need the gluon propagator 

5-'(Cvr)=E2(l-(-lf^). (5.35) 
p 

Altogether, after integrating over gluon momenta near k = Ctt to remove the delta-function 
in Eq. p.32| , we find the following non-zero matrix elements 



((7,c®Cpc))£ = 4;^,¥eB{Cn)\{-lf^-{-lf^]^=W^^,^^, (5.37) 



where fi, p and a can take any values, but are all different, and C is any non-zero hypercube 
vector satisfying = 0. The remaining one- loop contributions to p'"* can be read off from 
this equation. There are no continuum contributions to these coefficients, so pij = — p'j*. 

These results show the complications due to the use of staggered fermions. The mixing 
is between bilinears all of which have F ^ I; the flavor is transmitted by the highly off- 
shell gluon. Furthermore, the mixing of Eq. |5.37| is flavor off-diagonal, but is nevertheless 



consistent with the hypercube symmetry group BTI], as we have checked explicitly. 



5.5 Applications 

The continuum octet operators which are needed to calculate kaon decay amplitudes are 

Ql = Sa7M(l + 75)Mf) M67p(l+75)c?a , 

Q2 = SaJ^,{l + j5)Ua M67m(1 + 75)4 , 

Q3 = Sa7^(l+75)rf„ E ^67^(1+75)^6 , 

Q4 = Sa7;,(l + 75)4 E Qblf^{^+l5)qa , (5.38) 
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where we follow the notation of Refs. [25, 26|, except that we use (I + 75) for the the left- 



handed projector. The sum over q runs over the / active flavors: u, d, s and possibly c. For 
brevity we have left out the electromagnetic penguin operators, Qt-io- It is straightforward 
to extend the discussion to them. 
References \\ 
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have calculated the two-loop anomalous dimension matrices for these 
operators in the NDR scheme. (They also use the HV scheme which we do not consider 
here.) We wish to find the lattice operators which match onto Qi-q to one-loop accuracy. 
We consider here only the contributions to the matching from penguin diagrams; those from 
X diagrams can be obtained from sect. ^ or, in the case of gauge- invariant operators, from 
Ref. 



19 



The first step is entirely in the continuum. It turns out to be more convenient to change 
from the basis of operators given above to that in which Qi, Q2, Q3 and Q5 are replaced by 
their Fierz transforms 



Ql = Sa7^(l+75)(ia Mb7^(l +75)^6 , 
Q2 = 5^7^(1+75)4 Mfe7/.(l+75)Ma , 

Qs = J2 Sal,,{l+75)qb g67M(l+75)4 , 
q 

Q5 = -2^ Sa(l-75)gb gfe(l + 75)4 , 



(5.39) 



respectively. The two bases are identical in four dimensions, but differ in n-dimensions. In 
particular, the finite parts of their penguin matrix elements differ. The matching coefficients 
between these two bases are simple to obtain by generalizing the work of Ref. , and we 
give them below. 

When we transcribe these operators onto the lattice, we do so at the level of contractions 
1^, [1^. For example, the penguin contraction of Q2 in a i^' — > tttt matrix element is obtained 



using the penguin contraction of the lattice operator 



Nf 



iXs{l^^5 ® ^5)Xd Xuilf^ ® I)Xu)i + (Xs(7m ® ^5)Xd X«(7m5 ® I)Xu) 



(5.40) 



The 1/Nf cancels the extra lattice flavors in the loop. In all cases of practical interest the 
flavors of the bilinears are either / or ,^5 . This means that the mixings of Eqs. ^.36| and 
5.37| are not relevant, because the flavor of the bilinears in these equations is never I or ^5. 
The only relevant mixing is that of Eq. |5.24| , which affects just the bilinears (7^ /). Thus 



only the flrst term in Q2 receives corrections. 

We can now explain why we use Q2 rather than Q2. When we transcribe the former onto 
the lattice, as in Eq. |5.40| , the penguin loop attaches to a single bilinear, and the spin-flavor 
indices are contracted in two loops. On the other hand, had we transcribed Q2 onto the 
lattice, the spin-flavor indices would have been contracted into a single loop. While there is 
no fundamental reason to prefer one choice over the other, numerical simulations have so far 
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used the two spin-flavor loop form 0. Thus we have calculated the corrections in sect, ^for 
this form. Given this, we must match onto Q2 rather than Q2- A similar argument applies 
for Qi, except that the penguin contraction vanishes at one-loop because of the color factor, 
so it turns out to make no difference whether we choose Qi or Qi- 

For the operators Qs-e, the situation is more complicated. Each operator has two types 
of penguin contraction: either q is contracted with g, or the penguin loop is formed by 
contracting q with s or g with d. For each operator, one type of contraction has a vanishing 
color factor.^ To put the remaining contraction in two spin-flavor loop form, we must 
match onto Qs, Q4, Q5 and Qq. The lattice transcriptions of the penguin contractions of 
these operators are analogous to Eq. ^.4C , and we do not give explicit forms. 



We now collect the penguin contributions to the matching coefficients. For the negative 
parity parts of the operators, the general form is 

Q,(NDR) = + P[Q,]-|!^Q^-* , (5.41) 



Q'p = ll^E [fc(7/.5®e5)Xd X,(7m®^)x,)/- |fe(7M5®e5)Xd Xg(7M®^)x,; 



For the positive parity parts of the operators, the same result applies with (7^5 ® ^5) replaced 
with (7^ ® J) in QpK The calculation done in the previous two subsections applies directly 
to the operator Q2. For the other operators we need only count the number of continuum 
flavors that can run round the penguin loop. Including the matching coefficient between the 
continuum bases, we find 

P[gi] = , P[Q2] = \ + Ip, Pm = 2P[Q2] , Pm = PlQe] = -Jp , Pm = O , (5.42) 



where p = p^i^i- One check on Eqs. |5.30| and ^.41| is that we can extract the contribution of 



penguin diagrams to the one-loop anomalous dimension matrix from the coefficient of ln(/i). 



Our results agree with those of, for example, Ref. |^ . 

How large are the matching corrections? First we note that Qp = —Q3/S+Q4—Q5/3+Q5. 
Thus the operator produced by mixing is of the same "size" as the original operators. To 
estimate the size of the coefficients, we proceed as above and set fi = n/a and g^/lGrr^ ~ 1/90. 
The largest coefficient of Qp is {g'^/16'7i'^)fp/{2Nf) ^ —1/20. Thus the penguin contributions 
to matching are small. 



There is a subtlety in the application of Eq. |5.41| to the quenched approximation. Recall 



that the initial mixing is with a bilinear containing gluon fields, which is then rewritten as 
a four-fermion operator using the equations of motion. The second step is altered in the 
quenched approximation. It is simple to show that, when taking matrix elements of the 
operator produced by mixing one should exclude contractions between the Xq and Xg fields. 

5.6 Beyond one loop 

At higher orders in perturbation theory the continuum operators Qi-e will match onto a 
linear combination of a much larger number of d' = 6 lattice operators. We wish to point 



^^We stress that in an actual numerical simulation one must include both types of contractions — only at 
one-loop order does one of them vanishes. 
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out some general features of this mixing. Some of the results stated here rely on the analysis 
of sect. ^ 

Mixing produces both four-fermion and bilinear d' = 6 operators. We are only interested 
in those operators whose flavors match those of the external states, i.e. those having flavor 
(^i? = 1 or 75. One- loop mixing produces essentially all the four-fermion operators satisfying 
this criterion: higher loops hold no surprises. This is not true, however, for the bilinears. In 
particular, the lattice transcriptions of the operators 

^57^(1+75)/^/^/^^ and ^sB^B^B^7^(l+75)rf, (5.43) 

can appear. These are allowed by the lattice symmetries, though not by the continuum 
symmetries, and cannot be simplified using the equations of motion. Normally such operators 
are suppressed by powers of a, but this is not true here because they are corrections to 
d' = 4 operators which themselves are multiplied by a~^. The operators of Eq. |5.43| are not 



produced by one loop graphs of Figs. ^ and |^, because these do not depend on the momenta 
of the s and d quarks. Presumably they appear in two loop graphs such as that shown in 
Fig. 1(e). 

Because they appear only in two-loop graphs, these operators are likely to be unimportant 
numerically. We mention them mainly as a curiosity. We note, however, that they also occur 
with Wilson fermions, yet have so far been overlooked, while other operators appearing at 
two- loops have been considered |Tl| . The difference of the two operators in Eq. |5.43| is even 



under CPS symmetry (see Ref. and the sect. H), and so does contribute to both even 
and odd parity matrix elements. Thus the tricks suggested in Ref. 110] do not remove this 
operator. 

We mention also two other d' = 6 operators which appear at two (or higher) loops 

^9^[s7^F^^(l+75)rf] and ^ 9^ [57^75^/.:. (1+75 )c?] ■ (5.44) 

These are allowed by the continuum symmetries, and by CPS, but are usually dropped 
because, being total derivatives, they do not contribute to physical matrix elements. On the 
lattice, however, the standard approaches to calculating K ^ rni amplitudes use indirect 
methods based on amplitudes in which the weak Hamiltonian inserts momenta. For example, 
one method requires the K tctt amplitude at threshold for equal mass quarks, so that 
energy is inserted by the operator [|l^. Thus the operators of Eq. p.44| can contribute. 



5.7 Landau gauge operators 

We close this section with a brief discussion of the differences in the results for Landau-gauge 
and gauge-invariant operators. The only difference in the one-loop calculation is that Fig. 
0(b) is absent for Landau gauge operators, leading to a violation of the Ward identity, Eq. 
|5.17| . The contribution of this diagram to P^lk; B, A) is (for ma = 0) 

/ Si<P)BAiA-B),e''^<^-^^ = -An'J^i;;^!)^^ , (5.45) 

J (h 
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which is independent of k. By enforcing gauge invariance the Ward identity restricts the 
mixing of the bihnears. Eq. |5.45| shows that, without this restriction, there is unwanted 
mixing of the Landau gauge bihnear with the lower dimensional operator A^. Because of 
this, a continuum operator such as Q2 will mix with the lattice transcription of 37^(1 + 
75)A^(i. As mentioned at the beginning of this section, such mixing cannot be corrected 
using perturbation theory; it must be eliminated by a non-perturbative constraint. One 
must also worry about mixing with operators containing two or three factors of A^. Thus, 
Landau-gauge operators require several non-perturbative subtractions, which makes them 
unattractive for use in numerical calculations of penguin diagrams. 

Barring this caveat, the mixing of unsmeared Landau gauge operators with the penguin 
operator is identical to that for gauge-invariant operators. This follows since the equations 
of motion for the external fermion line annihilate any contribution from the longitudinal part 
of the gluon propagator. 

For smeared Landau gauge operators additional work is required, but, given the problems 
with gauge noninvariant operators, we have not done the calculation. 



6 LOWER DIMENSION OPERATORS 

We close the paper by considering the mixing of SU(3) octet operators with lower dimension 
operators. We begin by summarizing the results of Ref. p|] for continuum operators. The 
diagrams of Figs. and ^ lead to mixing of the d = 6 operators listed in Eq. ^.38| with 
quark bilinears of lower dimension, d' < 6. The original operators transform in the (8, 1) 
representation of SU{3)l x SU{3)ji, are Lorentz (pseudo)scalars, and have positive parity 
under CPS (CP followed by s 0? interchange). The only lower dimension bilinear with 
these properties is 

0,ub = S7m(1+75) Kd, 5= (/^ - 5^) . (6.1) 

In the continuum, one does not need to worry about this operator because, after using the 
equations of motion, it can be absorbed into the mass matrix by a chiral rotation of the 
fields. On the lattice, however, one must subtract Osub by hand This cannot be done 



using perturbation theory because the coefficient is proportional to 0. One must use 
a non-perturbative method such as that suggested in Ref. This uses chiral perturbation 
theory, and relies on the fact that the positive and negative parity parts of Osub are related. 
The remaining issue is whether lattice artifacts invalidate the subtraction procedure, for 
example by causing mixing with additional lower dimension operators. 

For Wilson fermions the subtraction procedure is, indeed, invalidated by the explicit 
breaking of chiral symmetry in the action. For example, there is mixing with the d' = 
3 bilinears sd and S'j^d, with no relationship between the two coefficients 0. Staggered 
fermions, on the other hand, do have a remnant chiral symmetry on the lattice, and it 
turns out that this is sufficient to allow the method of Ref. @] to be employed. Part of 
the argument has been given in Refs. [|T2], |T3|, |^, the remainder is given here. The 



argument only applies, however, if the original four- fermion operators are gauge-invariant, 
and we assume this henceforth. 

One point has been glossed over in this discussion. We have assumed that the operators 
produced by mixing have definite transformation properties under chiral SU{3), as is true 
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in perturbation theory. We must include non-perturbative effects, however, since the mixing 
is with operators of lower dimension. In particular, chiral symmetry itself is broken by 
the vacuum. Nevertheless, the assumption remains valid because the chiral transformation 
properties of the operators imply that their matrix elements satisfy Ward Identities, and 
these are not affected by the breaking of chiral symmetry. 

There are two questions that must be addressed with staggered fermions. First, if we 
transcribe the (8, 1) operators onto the lattice, using the rules explained in Refs. [|1^, |^, [l^, 
with which lower dimension operators can they mix? Second, are the positive and negative 
parity parts of the operators that one obtains related in such a way that one may use the 
subtraction method of Ref. H]? The answer to the first question, as we show below, is that 



the only relevant mixing is with the lattice transcription of Osub- Given this, Refs. |T^ , [T^ 
have answered the second question affirmatively. Thus in the following we discuss the first 
question alone. We consider only the positive parity parts of operators, i.e. those which 
contribute to the K it matrix elements. The extension to the negative parity parts is 
obtained by multiplying the operators listed below by (75 (g) ^5). 

Mixing is restricted by the properties of the operators needed to transcribe the penguin 
contractions of Qi-e onto the lattice. To illustrate these properties we show two examples. 
The lattice transcription for the positive parity part of Qi is 

NfQi' = {xsK^Xd Xuil^®I)Xu) + {Xsilf.5 ® I)Xd X„(7/.5 ® I)Xu) (6.2) 

Only the penguin contraction of this operator is to be considered, i.e. that in which Xu ^i^d 
Xu are contracted. For the operator the lattice transcription involves several operators. 
That giving the largest contribution is 

N^Q't = -2{xs'(T^Xd[xs{T^Xs + XaW^Xd])^ 

+ 2 (x.(75®e5)Xd [x.(75®6)Xs + XdVf^^^Xd] ) ^ . (6.3) 

Different contractions of the two parts of the operator are to be kept: for the first line, only 
those in which the quark fields of the second bilinear are contracted with each other; for the 
second line, only those where each bilinear is contracted with one of the external mesons. 
The second contraction is not a penguin diagram, but it is crucial that it be included in 
order to satisfy the Ward Identity described below. The need for this extra term is related 
to the fact that all the bilinears in Eq. |6.3| have even distance, in contrast to the odd distance 
bilinears in Eq. p.2| . 

These operators, and all others which arise, have the following properties: 

1. They are singlets under the group of reflections and rotations which map the hypercube 



into itself, the hypercubic group |Q. Operators produced by mixing must also be 
singlets. 

2. They are flavor octets. Thus the bilinears produced by mixing must have the flavor 
form Xs{ls®^F)Xd- 

3. They have positive parity under CqS, where Cq is the lattice charge conjugation symme- 
try (see, e.g. Ref. |0), and 5* is the s ^ d interchange symmetry. The transformation 
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rules for CqS are 



[Uplab ^ [Ul]ba , [Fpu]ab ^ -[Fpu]ba , "^s ^ "^d, 

fc]a(75®eF)M5 ^ (-)^[xj6(4®ef)Ma , (6.4) 

where a, b are color indices, and F^j^^, is a lattice transcription of the gluon field strength. 
The lower dimension operators must therefore have positive CqS* parity. This is possible 
either if the bilinear itself has positive CqS parity, or if it has negative CqS* parity and 
is multiplied by — m^. 

4. They satisfy a Ward Identity relating their K^tc^ correlator {Cxn) to that between the 
strange scalar k = sd and the vacuum (C^): 

\/2{ms + md) CxnitK, U) = -C,,{U) . (6.5) 

(This is equivalent to Eq. 8.10 of Ref. |]12|].) This relation holds provided that one uses 
the operators (75 ® ^5) and (/ ® /) to create the pseudoscalars and scalar, respectively, 
and provided that m„ = m^.Q It is the lattice analog of the continuum Ward Identity 
which follows from current algebra. It is straightforward to show that Eq. ^3] can be 
satisfied by lower dimension operators in two ways: (a) if the bilinear itself has positive 
CqS parity, then it must have odd distance; (b) if the bilinear itself has negative CqS 
parity, then it must have even distance. These selection rules play the role of parity in 
the continuum CPS symmetry. 

Verstegen has classified staggered fermion bilinears under the hypercubic group pT|.p^ 
Using his results we find that the following lower dimension bilinears have all the symmetry 
properties listed above and thus can be produced by mixing 



dim3 : S3 = J] X.(75 ® ^p5)Xd 



dim4 : Bia = Yl Xsilf^ ® I) Dp Xd 



dimS : iSg^ = ^ Xshp<T^Cp5) Fpi/Xd 



pupcr 



Xs{lpa®^p^) 



pupa 



135c = J2XsDp (75 ® ^p5) Dp Xd 



IS^d = J2XsDp (75 ® ^p5) D^ Xd 



^•^For the even distance operators such as that in Eq. there is a subtlety concerning the allowed 

contractions of the operator in Ck which is discussed in Ref. jlj]. Here we simply assume that Eq. [b!^ 
applies. 

^"^In applying his results one must note that his flavor matrix corresponds to our ^^5. 
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V Xd 

pupcr 



pupa 



Xd 



= J2 ^^^^ypcr^p[Xsi'yp<T ® ^p5) Du Xd] 



fii/pa 



f Xd\ 



pupa 



B,^ = ^'J2[xAl5^U)Xd] 



B5k = E W'x^M'9z.[Xs(7pa 



pupa 



F^y is any lattice version of the field strength with appropriate transformation properties 
under the hypercubic group. The covariant derivatives will be defined precisely below. 
We have included total derivatives in the list because they can contribute if one calculates 
K ^ IT matrix elements with non-degenerate quarks. All the operators given above can 
be multiplied by functions of nis and rrid that have positive CqS parity. It turns out that 
no lower dimension operators with explicit factors of — rrid (or any function of rris and 
rud having negative CqS parity) satisfy all the properties listed above. This is why all the 
operators in the list have odd distance. 

The original four-fermion operators (e.g. Eqs. and |6.3| ) are staggered flavor singlets.Q 
Thus, in the continuum limit, only bilinears with F = I can be produced by mixing. In fact, 
of the fourteen lower dimension operators, only B^a has F = L It is the lattice transcription 
of Osub, and is the operator we expect to be produced. All the other operators are staggered 
flavor non-singlets, and are allowed because elements of the hypercubic group contain flavor 
transformations. What we must determine is the conditions for which these operators, which 
are artifacts of the lattice regularization, do not contribute to matrix elements. 

The first condition is that we take matrix elements between states having the same 
staggered flavor, e.g. between two pseudo-Goldstone pions each having flavor ^5. At first 
sight this solves the problem because the matrix elements of staggered flavor non-singlet 
operators should vanish, leaving only B4a- This is only true, however, up to corrections 
of 0{pextO'), where pext is an external physical momentum. In effect, the matrix element 
picks out higher dimension parts of the original operator, and these parts can have different 
staggered flavors. What one must do is take out sufficient factors of momenta until one 
obtains a flavor singlet. The flavor singlet components of all the non-singlet operators listed 
above turn out to be of dimension 6 or higher. 

To explain how this works we need some results concerning the full staggered fermion 
symmetry group including translations [Q. In particular, we need the representations of 
the zero four- momentum subgroup, Qq, which are summarized in Ref. |]T2|. To discuss these 
representations, it is convenient to use a different definition of bilinears. Our choice up to 



15 



We use the expression "staggered flavor" to distinguish from the continuum flavor SU(3) symmetry. 
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now, Ospiy), (defined in Eqs. and |2 .121) has all fields lying in a single hypercube, labeled 
by y. This definition has the disadvantage that Y^y^sriy) is not a representation of ^o, 
except for operators with A = 0. 

Instead, following Ref. we proceed as follows. Each term in Osriy) is a product 

xiy + ^)x{y + B) multiplied by a phase, where B =2 A + S — F , the subscript indicating 
mod-2 arithmetic. We first think of this term as if both x aiid x started at site y + A, with 
X then shifted to site y + B. If (S* — F)^ 7^ 0, then in the new operators we perform this 
shift symmetrically 

S^x{y + A) = \[x{y + A + li) + x{y + A - ii)] , (6.6) 

rather than picking just one of these two terms. We do this symmetric shift in all directions 
for which 5* — F is non-vanishing, resulting, in general, in an operator contained in a 4*^ 
hypercube. To complete the definition we take the average of this enlarged operator with 
that obtained by shifting x instead of x- This last step gives an operator which transforms 
into itself under Cq, with the sign given by Eq. The resulting operator we call Ogp'^iy), 
because Y^y^tWiv) is part of a true representation of Qq. In particular, it transforms into 

itself under translation in the /x'th direction, with phase (— )^'', where F^ =2 X^i/^^-^i/- 

We also need additional "derivative" operators defined as for O^jf except that one or 
more of the symmetric shifts is replaced by the antisymmetric shift 

D^xiy + A) = llxiy + A + fx)-x{y + A-fi)]. (6.7) 

This corresponds to taking a derivative of the quark field, and we write the resulting bilinear 
ci-s xils ® ^f)D^x- Clearly {S — F)^ must be non-zero for there to be the possibility of an 
insertion of D^. Unlike for Og'p'^, we distinguish explicitly between the derivatives acting on 
the quark and antiquark fields. These must be combined to obtain eigenstates of CqS, as in 
B^a- Two properties of derivative operators are important. First, when summed over all y, 
they give true representations of Qq, with the translation phase being (— )'^''. Second, they 
transform under the hypercubic subgroup in the same way as x(7aj55 ® ^Ai5F)x- This follows 
because a hypercube operator Osf can be written as a sum of Og'p'^, which transforms exactly 
as Osf under the hypercubic group, and a number of derivative operators. For example. 



X(75 ® U)X = [X(75 ® - ®I)D,X- (6.8) 

This shows the decomposition of the original hypercube operators into true representations 
of Qo- 

We now return to the list of lower dimension operators Ei. First we complete their 
definition. Some of them contain covariant derivatives in directions for which — F = 0, 
e.g. B^c- In these cases the covariant derivative is defined by 

D^x{y + A) = \ [xU + A + 2P)- xU + A- 2/2)] . 

Similarly the total derivatives are defined by shifting two lattice units, which avoids mixing 
with staggered fiavor symmetries. Next we note that, using results such as Eq. |6.8| , we can 
rewrite all the operators so that they become true representations of Qo when summed over 
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y. In fact, the list already includes all the derivative operators of dimension less than 6 that 
can appear. For example, turning Bz into a true representation using Eq. |678| , we find the 
derivative operator B^a- Thus, with no loss of generality, we can take all the bilinears in the 
list to already be true representations. 

We can now make more precise the statement that the matrix elements of the flavor 
non-singlet operators vanish up to momentum corrections. Since F = iff F = 0, all such 
operators have a non-trivial translation phase in at least one direction. Thus their matrix 
elements vanish aX p = 0, where p is the total momentum inserted by the operator. In 
order to obtain, say, a K it matrix element aX p = 0, the K and vr must both have the 
same staggered flavor (since otherwise p ~ Ctt for some non-zero hypercube vector C), and 
the K and vr must be degenerate (for one can then create the states at infinite positive and 
negative time, respectively, so that the operator can be averaged over all positions). For such 
a matrix element only B^a contributes and the argument is complete. In a general lattice 
matrix element, however, there will be a small physical momentum inserted by the operator, 
e.g. if the K and tt are not degenerate. In this case the flavor non-singlet operators can have 
non- vanishing matrix elements. 

To evaluate these contributions we use two identities concerning matrix elements when 
(-)^f' = -1. The first is vahd when {S - F)^ ^ 0, 

((x(75®eF)x)*™^(p)) = tan(p,,/2)2(lx(7M5S®W) K x) , (6.9) 
the second for {S — F)^ = 

((x(75®eF)x)*™'^(p)) =aan(p^/2)(x(7M55®UF)x) • (6.10) 

There is no implied sum over /i in either equation. The operators are defined on a single 
hypercube, and the matrix element is between external states which insert momentum p. 
For small p, the factors of itan(p^/2) can be replaced by acting on the operators. 

These equations remain valid even if there are covariant derivatives or field strengths in the 
operators. Since (— )'^^^m = +1, the operators on the right-hand-sides have trivial translation 
phases in the /i'th direction. By using the equations sequentially in all directions in which 
the translation phase is non-trivial, one can find the flavor singlet component in a given 
bilinear. 

We now apply these equations to the bilinears Bi. Each application raises the dimension 
of the operator by either 1 or 3, depending on whether {S — F)^ = or 1. Recall that all 
the bilinears Bi are now taken to be true representations of Qq. We work our way up from 
the (i' = 3 operators. 

d'=3 There is only one such operator: B^. It can be replaced by \^'^Xsi^^^ ® ^) -^m Xd, using 
Eq. |6.9| . This is a dimension 6 operator which, furthermore, vanishes by the equations 
of motion if the quark masses are set to zero. 

d'=4 The only flavor non-singlet dimension 4 operator can be written 

= E d,B^ ' = E xJa^Mi:^Xd . (6.11) 
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One application of each of Eqs. |6]^ and |6.10| shows that acts as a dimension 7 



operator in matrix elements, so that 841, acts as a dimension 8 operator and can be 
ignored. 

d'=5 Dimension 5 operators are converted into flavor-singlets of dimension 6 or higher using 
Eqs. |6.9| and |6.10| . All dimension 6 operators that occur will appear in the mixing 
calculations at some order. In sect. |S| we calculated this mixing to one-loop, and 
found that the only staggered flavor singlet operators that appear can be rewritten as 
four-fermion operators. 

We conclude that all the lower dimension operators except B^a have vanishing matrix ele- 
ments even if non-zero physical momentum is inserted. 

We close by noting that the term proportional to /„ in Eq. p.l8| is an explicit example 



of mixing with lower dimension operators. At one loop this does not affect the lattice 
transcriptions of Qi-q, but when inserted in a two loop diagram, it gives rise to mixing with 
the operator B^a- This operator has (— )^'' = —1, and {S — F)p 7^ 0. Using Eq. |6l9| , one 
finds that it acts in matrix elements as though it has dimension 8, and thus can be ignored. 
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A FIERZ TRANSFORMATIONS 

In this appendix we give the Fierz transformation rules for the 35 diagonal scalar opera- 
tors. We do not include the sign due to Fermi statistics. Since we always do two Fierz 
transformations, any such sign would cancel. 

The first 25 operators are of the form exemplified by [V" x T], in which the spin and 
flavor matrices are independent. These form a closed set under Fierz transformation, with 
the spin and flavor part of the operators being transformed independently. The rules can 
thus be obtained from those for the spin alone: 

V-A ^ A-V 

V + A ^ 2{S-P) 

S + P ^ i(5 + P + T) 

T ^ i(35 + 3P-T) . 



Thus, for example, the transformation needed in the example of Section 3.3 is 



[V X S] ^ i [{2S -V + A-2P)x{S + V + T + A + P)] . 
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The Fierz transformation rules for the remaining 10 operators are 



[{y + A)^ x{y + A)^] ^ [{S -P)x{S- P)] + [T_ X T_] 

2[{V + A),xiV-A),] ^ -[(S-P)xiV-A)]- [T,xiV + A),] 

[(V~A),x{V~A),] ^ [(V-A),x{V~A),] 

[{V~A),xT,] ^ [T,x{V-A),] 

2[{V + A)^xT^] ^ [(V-A)x{V + A)]~A[iV-A)^x{V + A)^] 

4 [r_ X r_] ^ - [{V + A)x{V + A)]+A [{V + A),, x{V + A),,] 

4 [T+ X T+] ^ 4 [{S + P) X {S + P)] - [{S + P + T) X {S + P + T)] 
+ 4[T+xT+] . 



B Corrections from Xa diagrams 

In this appendix we use the one-loop matching coefficients for bilinears to calculate those 
for the diagonal scalar four-fermion operators. We consider only the contributions from Xa 
diagrams, since, as discussed in the text, we can bring all other one-loop calculations into the 
form of an Xa diagram. We leave out color factors and an overall g'^/lGn'^, and set fi = n/a, 
so that the anomalous dimension term is absent. Thus, in the notation of Eq. |3.11| , we are 
calculating the matrix Ai"^. 

The matching coefficients for bilinears (defined in Eq. |3.8| ) are given in PS. They depend 
on both the spin and flavor of the bilinear, and this is codifled in the labels given to the 
Cj. For example, the correction for the bilinear {•y^^ (g) /) is labeled cas- Since cas does not 
depend on the index u, the correction for the operator [A x S] is diagonal. In the notation 
of Eq. |1| 

M^i [A X S]) = 2cAs [AxS] . (B.l) 

The factor of two comes from the fact that both bilinears are corrected. 

An important feature of the one-loop corrections for bilinears is that they commute with 
multiplication by (75 (8)^5). Thus if one simultaneously changes the spin and flavor of the 
operators according to S^P,V^ A, T^T, the correction coefficients remain unchanged. 
This means, for example, that cas = cvp- As in sect. ^ we define V to be the operation 
of multiplying both bilinears in the four-fermion operator by (75 ® ^5) from the left. Due to 
the axial symmetry, this operation commutes with multiplication by the correction matrix 
For example, applying V to Eq. |B.1| we find 



M\ [V X P]) = 2cAs [VxP] = 2cvp [VxP] . (B.2) 

The simple diagonal form of the correction in Eqs. [B.l| and [B.2| applies for those operators 
whose flavor or spin is either S or P. In this case, the corrections do not break the flavor or 
spin symmetries. For other operators, e.g. [l^ x "K], these symmetries are broken. This is 
because the bilinears in these operators consist of parts of different distances. The bilinears 
in [V X V] can have A = 2 (e.g. (7^ <S) ^u) with /i 7^ z/) or A = (e.g. (7^ <S) ^^)), and these 
are corrected with different coefficients {cvv2 and cyvo respectively). The net result is that. 
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in terms of the operator basis that we use, the corrections are 



M'^ilVxV]) = 2cvv2[VxV]+2{cvvo-cvv2)Ky<V^] (B.3) 

K X V,]) = 2cyyo K X ^m] • (B-4) 

Similarly, we find 

M''{[VxA]) = 2cvA2[V X A] + 2{cvA4 - cvA2)Ky< A^] (B.5) 

M'^i [V^ X A^]) = 2CVAA [V^ X A^] . (B.6) 

The remaining results for operators containing vector bilinears are 

M\[VxT]) = {cvTi + cvT^)[V xT] + {cvTi-cvTz)[V^^T^] (B.7) 
M\[V,xT^]) = {cvTi-cyTz)[VxT] + {cvTi + cvT^)KxT^]. (B.8) 

The results for operators containing axial bilinears can be obtained by applying V to the 
above results. 

Finally, we consider tensor operators. The results for [T x V] and [T^ x V"^] are as in 
Eqs. |B.7| and p.8| , with V and T interchanged. The relevant coefficients are thus ctvi and 
CTV2,- Results for [T x A\ can be obtained using V. The remaining results are 

M\[TxT]) = {ctt^-ctta)[T-xTJ\ + {ctto + ctta-2ctt2)[T+xT+] 

+ 2cTT2 [T X T] (B.9) 
M\[T_xT_]) = (c^^o + crT4) [T_ X T_] + (ctto - CTT4) [T+ X T+] (B.IO) 

M\ [r+ X T+]) = (CTTO - CtTa) [T- X r_] + {CTTO + ctta) [T+ X T+] . (B.ll) 
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